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The  Fleet’s  interest  la  the  Arctic  environment  is  constantly  increasing,  and 
the  need  to  be  able  to  understand  and  forecast  the  arctic  environmental 
phenomena  has  become  an  taportani  issue  of  ocean  research.  This  report  deals 
with  a  one-dimensional  coupled  ice/ocean  model  that  embodies  the  relevant 
Arctic  thermodynamic  and  dynamic  phenomena.  This  work  is  a  first  step 

toward  a  three-dimensional  ice/ocean  model.  ) 
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Executive  summary 


A  coupled  one-dimensional  ice/ocean  model  is  developed  in  differen  ial 
form.  The  ice/snow  system  is  represented  by  the  simplified  thermodynamic 
model  of  Semtner  and  a  dynamic  approximation  that  neglects  the  internal 
stresses.  The  ocean  is  represented  by  the  Mellor-Yamada  Level-  2  turbulence 
mixed-layer  model. 

The  thermodynamic  coupling  considers  a  moving  ice/ocean  interface  and 
a  salinity  flux  generated  by  the  freezing  or  melting  of  ice.  The  dynamic  coup  ing 
occurs  via  the  turbulent  stress  that  exists  in  the  mixed  layer  beneath  the  ice. 

Two  test  cases  are  used  for  model  validation  and  scientific  studies.  One 
is  the  standard  climatological  test  used  by  Semtner  and  others.  The  other  test 
case  is  with  the  A1DJEX  data. 

The  ice/ocean  model  is  compared  to  Semtner’s  ice  model  to  determine  the 
effect  of  a  variable-depth  mixed  layer  as  opposed  to  an  isothermal,  fixed-depth 
mixed  layer.  The  sensitivity  of  the  ice/ocean  model  to  changes  in  certain 
parameters  and  forcing  is  also  examined.  Finally,  the  ability  of  the  model  to 
simulate  Arctic  conditions  is  shown  by  comparison  of  model  results  to  data 
obtained  during  AIDJEX. 
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1.  Introduction 

Our  work  was  motivated  by  the  need  to  further 
develop  the  coupling  of  Arctic  ice  and  ocean.  The  pres¬ 
ent  status  of  the  coupling  ranges  from  an  ocean 
described  by  prescribed  heat  fluxes  and  currents  to  a 
■oarse  coupling  of  ice  and  ocean.  We  will  introduce 
a  detailed  polar  mixed-layer  representation  and  couple 
the  Arctic  ice  and  ocean  through  it.  This  will  enable 
the  ice  to  influence  the  ocean  and  the  ocean  to  influence 
the  ice.  We  will  develop  and  validate  the  representa¬ 
tion  of  the  associated  physical  processes  and  coupling 
mechanics. 

Most  of  the  previous  ice  models  treated  the  ocean 
in  a  passive  or  a  semipassive  sense.  The  passive 
approach  consisted  of  a  prescription  of  climatolog¬ 
ical  ocean  currents  and  heat  fluxes.  The  semipassive 
approach  included  a  heat  budget  calculation  with  a  con¬ 
stant  depth  mixed-layer  prescription  (Hibler,  1980). 
The  climatologies  ranged  from  mean  annual  to 
seasonal  time  scales.  Examples  are  the  work  of  Coach¬ 
man  and  Aagaard  (1974)  for  annual  ocean  currents, 
or  the  diagnostic  use  of  the  Hibler-Bryan  (1984)  results 
for  providing  seasonally  varying  ocean  currents  and 
heat  fluxes.  In  ice  forecasting,  the  use  of  the  seasonally 
varying  ocean  currents  and  heat  fluxes  has  been  shown 
to  yield  a  more  realistic  ice  edge  (R.H.  Preller,  per¬ 
sonal  communication,  1986). 

The  coupling  of  ice  and  ocean  has  so  far  been 
approached  at  various  levels  of  sophistication.  The 
thermodynamic  part  usually  contains  Semtner’s  (1976) 
snow  and  ice  model.  The  simplest  thermodynamic 
coupling  to  the  ocean  involves  the  prescription  of  a 
constant  depth  mixed  layer  and  a  heat  budget  calcula¬ 
tion  (Parkinson  and  Washington,  1979;  Hibler  and 
Bryan,  1984). 

More  sophisticated  thermodynamic  models  of  the 
snow,  ice,  and  water  system  involved  representations 
of  the  polar  mixed-layer  region  at  various  levels  of 
sophistication.  Pollard  et  al.  (1983)  have  experimented 
with  a  slab  model  superposed  on  a  prescribed  thermo- 
cline.  The  model  was  used  to  simulate  the  thermody¬ 
namic  cycle  of  a  snow,  ice,  and  water  system. 

The  simplest  dynamic  coupling  of  a  snow,  ice,  and 
water  system  was  accomplished  through  the  use  of  a 
drag  law  between  ice  and  ocean  (Parkinson  and 
Washington,  1979;  Hibler  and  Bryan,  1984).  A  more 


sophisticated  approach  incorporates  boundary  lajer 
formalism  and  uses  an  iteration  technique  to  calculate 
the  Reynolds  stress  between  the  ice  and  water  (Overland 
et  al.,  1984). 

Recently,  Mellor  et  al.  (1986)  developed  a  time- 
dependent,  second-moment,  turbulence  closure 
representation  of  the  coupled  ice  and  water  system. 
The  model  has  a  detailed  boundary  layer  representa¬ 
tion  and  considers  freezing  and  melting  phase  transi¬ 
tions,  and  represents  the  associated  unstable  and  stable 
states. 

Our  approach  considers  a  time-dependent,  horizon¬ 
tally  homogeneous  snow,  ice,  and  water  system.  The 
coupling  of  the  snow,  ice,  and  water  system  is 
developed  on  a  stretched  vertical  grid.  Second-order 
closure  is  used  for  the  representation  of  the  polar  mixed 
layer. 

The  thermodynamics  of  the  ice  and  snow  system  is 
represented  by  the  Semtner  0-layer  model.  The  ice 
dynamics  is  represented  by  momentum  equations  t  hat 
neglect  the  internal  ice  stresses.  The  coupling  of  the 
ice  and  polar  mixed  layer  is  developed  through  the 
interfacial  stresses,  heat  fluxes,  and  freezing/mel:ing 
processes  representation. 

Studies  of  the  climatological,  seasonal,  and  diurnal 
cycles  are  performed  with  the  coupled  ice/ocean  mcdel. 
The  previous  climatological  simulations  of  Semtner 
(1976)  are  repeated.  Feedback  due  to  ice/ocean  cou¬ 
pling  plays  a  role  in  these  simulations.  The  effecis  of 
feedback  can  be  seen  in  the  multiyear  equilibrium  cycle 
of  ice  thickness,  heat  budget,  surface  temperature,  and 
other  variables. 

The  cycle  of  open-water  appearance  changes  from 
Semtner’s  (1976)  previous  results.  Sensitivity  studies 
of  the  appearance  and  disappearance  of  open  water 
are  conducted  with  the  coupled  ice/ocean  model  by 
varying  various  parameters.  During  the  open-v.ater 
period,  the  water  is  heated  by  the  incoming  solar  radia¬ 
tion.  As  the  water  refreezes  a  hot  spot  moves 
downward  into  the  ocean  and  gradually  diffuses. 

Another  set  of  studies  is  performed  with  the 
AIDJEX  data.  The  data  spans  a  period  of  a  year  and 
contains  daily  measurements.  Model  simulations  of 
the  seasonal  cycle  are  performed  with  a  time  step 
that  resolves  the  diurnal  variations.  The  simulated 
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conditions  are  compared  with  the  AIDJEX  data.  The 
diurnal,  seasonal,  and  climatological  responses  are 
studied. 


17.5,  25,  32.5,  40,  50,  62.5,  75,  100,  125,  150,  2CO, 
300,  400,  and  500  m.  The  vertical  eddy  fluxes  are 
defined  midway  between  these  depths. 


2.  Description  of  the  model 

a.  Introduction 

This  paper  is  concerned  with  the  design  and  testing 
of  a  one-dimensional  numerical  ice/ocean  model.  In 
this  model,  Semtner’s  (1976)  0-layer  ice  model,  together 
with  simplified  ice  dynamics,  is  coupled  with  an  ocean 
mixed-layer  model  that  uses  a  turbulence  parameteriza¬ 
tion  scheme  based  on  Mellor-Yamada  Level-2  tur¬ 
bulence  closure  theory  (Mellor  and  Yamada,  1974). 

Semtner’s  ice  model  assumes  a  state  of  conductive 
equilibrium  for  heat  and  has  three  prognostic  variables: 
snow  thickness,  hs;  ice  thickness,  /»,,  and  surface 
temperature,  Ts.  His  model  uses  the  imbalance  in 
atmospheric  fluxes  at  the  snow/ice  surface,  and  the 
conductive  heat  flux  within  the  snow/ice  system  to 
compute  the  surface  temperature  of  the  system.  If  this 
temperature  rises  above  the  freezing  point,  the  surface 
temperature  is  set  at  the  freezing  point,  and  the  amount 
of  surface  melt  is  computed  from  the  heat  surplus.  At 
the  bottom  of  the  ice,  the  imbalance  between  a  pre¬ 
scribed,  constant  oceanic  heat  flux  and  the  conductive 
flux  within  the  snow/ice  system  is  used  to  determine 
the  accretion  or  ablation  of  ice.  Semtner  assumes  a 
30-m-deep,  vertically  isothermal  mixed  layer  and 
calculates  the  associated  heat  budget  during  open-water 
periods  in  order  to  simulate  seasonal  transition  between 
an  ice-covered  ocean  and  an  ice-free  ocean.  When  the 
mixed-layer  temperature  drops  to  or  below  the  freez¬ 
ing  point,  ice  reforms. 

A  simplified,  one-dimensional  ice  dynamics  model 
is  constructed.  The  model  assumes  a  horizontally 
homogeneous  layer  of  ice  and  neglects  the  internal  ice 
stresses.  This  is  the  free-drift  approximation  that  is 
valid  for  ice  drifting  in  the  ocean  and  not  experienc¬ 
ing  internal  stresses  and  pressure. 

The  vertically  isothermal  mixed  layer  in  Semtner’s 
model  is  replaced  by  a  dynamic,  thermodynamic  mixed 
layer.  This  polar  mixed  layer  is  represented  by  a  level-2 
turbulent  closure  model  (Mellor  and  Durbin,  1975; 
Clancy  and  Martin,  1981).  This  model  has  prognostic 
equations  for  temperature,  salinity,  and  velocity.  The 
model  is  in  differential  form  and  resolves  the  Reynolds 
and  thermal  stresses. 


c.  Initial  conditions  and  forcing 

Two  different  test  cases  are  used  for  model  valida¬ 
tion  and  scientific  studies.  One  is  the  standard  cent  al 
Arctic  climatological  test  case  used  by  Semtner  (19''6) 
and  others.  The  other  case  is  the  AIDJEX  experiment. 
The  AIDJEX  experiment  resolves  the  diurnal  time 
scale. 

For  both  test  cases,  the  mixed  layer  was  initialized 
with  temperature  and  salinity  profiles  taken  from  AID¬ 
JEX  data.  In  particular,  profiles  from  station  Blue  Fox, 
January  1,  1976,  were  chosen.  These  profiles  repre¬ 
sent  typical  winter  profiles  for  the  Beaufort  Sea,  where 
the  AIDJEX  data  was  taken  (Figs.  2  and  3).  For  :he 
climatological  test  case,  the  model  was  forced  with 
radiative  fluxes  from  Fletcher’s  monthly  mean 
climatology  (Table  1,  Semtner,  1976).  Semtnct’s 
snowfall  rate  and  average  monthly  winds  from  1983 
Fleet  Numerical  Oceanography  Center  (FNOC)  analy¬ 
sis  were  also  used.  Figures  4-10  show  these  fluxes  in  er- 
polated  to  each  time  step.  (The  AIDJEX  test  case  will 
be  described  later.) 

d.  Ocean  currents 

Ice  and  ocean  velocities  are  calculated  and  used  in 
estimates  of  mechanical  stirring  within  the  mixed  layer. 
A  constant  geostrophic  current  of  ug  =  -1.2  cm/sec 
and  vg  =  -0.65  cm/sec  is  assumed.  Horizontal 
advection  however,  is  not  included. 


e.  Ice  thermodynamics  model 

A  complete  description  of  Semtner’s  0-layer  ice 
model  is  contained  in  his  paper.  The  following  i  >  an 
outline  of  the  relevant  equations. 

The  conductive  heat  flux  through  the  snow/ice 
system  (assuming  heat  flux  and  z  are  pos  live 
downward)  for  snow-covered  ice  is 


-  ki(TB 
h.  + 

l 


Ts> 


(h  k . /k  ) 
v  s  1  s 


(1) 


and  for  snow-free  ice  is 


-  ki<TB 


V 


.r . 

l 


(2) 


b.  Grid 

The  vertical  grid  used  in  the  coupled  ice/ocean  model 
has  two  points  in  the  snow/ice  system  and  17  levels 
between  the  ocean  surface  and  500  m  depth  (Fig.  1). 
Temperature,  salinity,  density,  and  velocity  for  the 
water  column  are  defined  at  depth  of  2.5,  7.5,  12.5, 


where  k,  =  4.86E-3  cal/cm/sec/°C  and  k,  =  7.4E-4 
cal/cm/sec/°C  are  the  conductivities  of  ice  and  snow , 
respectively.  The  quantity  TB  is  the  temperature  at  the 
bottom  of  the  ice,  Ts  the  temperature  at  the  surface 
of  the  snow/ice  system,  h,  the  ice  thickness,  and  hs 
the  snow  thickness. 
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Table  1 .  Listing  of  parameters  and  assigned  values  used  in  Case  1  for  both 
the  CML  case  and  the  VML  case. 


Constant  Mixed  Layer 

Variable  Mixed  Layer 

H,: 

50  cm 

50  cm 

H.: 

30  cm 

30  cm 

F8: 

- 1 .929E-4  cal/cm2/sec 

- 1 .929E-4  cal/cm2/sec 

Ts: 

-2.0°C 

-2.0°C 

Tb: 

-2.0°C 

Varies  with  salinity 

o: 

1.385E-12  cal/cm2/K4/sec 

1.385E-12  cal/cm2/K4/sec 

a: 

0.66448 

0.66448 

<*«■ 

0.07 

0.07 

<v 

72  cal/cm3 

72  cal/cm3 

64  cal/cm3 

64  cal/cm3 

r ■ 

1.065 

1.065 

MLD: 

30  m 

Variable 

SFF: 

None 

3.341  E-5  cm/sec 

Time  step: 

8  hours 

8  hr/30  min 

Forcing: 

Heat  fluxes: 

Climatology 

Climatology 

Winds: 

None 

FNOC  average  monthly 

T&S  profile: 

None 

AIDJEX  data  (Beaufort  Sea) 

* 


> 
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The  change  in  the  surface  temperature  Ts  is  com¬ 
puted  from  a  heat  balance  between  the  atmospheric 
fluxes  and  the  conductive  flux  Fs  at  the  surface  of  the 
snow/ice  system.  The  surface  temperature  at  a  previous 
time  step  is  denoted  by  Tp  and  Ts  =  TP  +  ATS  is  the 
temperature  at  the  present  time  step.  Thus  ATS  can  be 
solved  from  the  heat  balance  by  making  the  above 
substitution  for  Ts.  The  result  is 

aTs  •  [F1  *  F,  *  FL  *  <1'«s>fr 

-  V*  ki<v 

*  (hski/ks»J 

'[‘Tp’  *  V<hi‘  <hski/ks»]'  (3) 

when  snow-covered,  and  by 

ATS  *  [F1  +  Ft  +  FL  +  (1_<xi)FR 

-  V  +  ki(TB-  V/hi] 

/[^?  -  kj/hjJ, 

(4) 


when  snow-free.  The  quantity  F ,  represents  latent 
heat  flux,  F,  sensible  heat  flux,  FL  incoming  longwave 
radiation,  a,  and  as  the  albedos  of  ice  and  snow, 
respectively,  FR  incoming  solar  radiation,  and  o  the 
Slefan-Boltzmann  constant.  Second-order  or  gret.ter 
terms  in  A Ts  are  ignored  in  the  previous  equations. 

The  accumulation  of  snow  is  determined  by  a  set 
snowfall  rate.  The  amount  of  snow  melt  is  computed 
from 

=  At  (FS  "  FA)/  V  (5) 

when  Ts  is  above  freezing. 

The  growth  or  decay  of  ice  at  the  bottom  is  deter¬ 
mined  from 

4hb  ■  «  <FB  -  FS)/  Ob-  (6, 

and  the  melting  of  ice  at  the  surface  when  snow-free 
and  Ts  is  above  freezing,  is  computed  from 

Ahf  -  At  (Fs  -  Fa)/  q..  (7) 

The  quantity  FA  represents  the  sum  of  the  atmos¬ 
pheric  fluxes,  i.e.,  the  first  five  terms  on  the  right  side 
of  Eq.  (3)  and  Eq.  (4),  FB  the  oceanic  heat  flux,  q  the 
volumetric  heat  of  fusion  for  the  snow,  qh  and  q  the 
volumetric  heats  of  fusion  at  the  bottom  and  the  top 
of  the  ice,  and  A t  the  time  step. 
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Ice  dynamics  model 

The  equations  defining  the  ice  drift  are 

u»  c .  _  v  H8 

W  S  "  _KH  3z  ’ 

(15) 

3u .  , .  .  x  T 

571  =  t(v.  -  v  )  +  -ax  +  -vx  , 

3t  1  g  m  ra  ’ 

(8) 

v'u'  =  -  K 

v  3z  ' 

(16) 

3V  •  r,  .  t  T 

t— 1  =  -r(u.  -  u)  +  -ay  +  -vy  , 
3t  v  1  g'  m  }  m  J  ' 

(9) 

v'v'  =  -K 

v  3z  ' 

1  17) 

where  u,  and  v,  are  the  x-  and  y-components  of  the  ice 
velocity,  and  ta>,  are  the  x-  and  y-components  of 
the  wind  stress,  twj  and  twv  are  the  x-  and  y-compo- 
nents  of  the  water  stress,  and  m  is  the  mass  of  the  ice. 


where  KH  and  KM  arc  eddy  diffusion  coefficients. 

For  a  more  detailed  description  of  the  ocean  model 
equations  see  Appendix  A. 


g.  Ocean  model 

The  ocean  model  consists  of  a  mixed  layer  and  the 
upper  ocean  down  to  500  m.  Physical  processes 
included  are  vertical  mixing,  radiation,  and  planetary 
rotation.  The  conservation  equations  for  temperature, 
salinity,  and  momentum  are 


It  =  h  (*^)+  h  <-v't,+  v  & 


i_ 

PoC 


3F 

3z 


SR 


(10) 


H  -  h  <-”s>  *  I?  v  S 


(ID 


=  fv  -  fv  +  |r  (-^r  +  V  f^v)  ,  (12) 
at  w  g  3z  w  az 

=  -fu  +  fu  +  ^_  (-v'v'  +  X)  ^v),(13) 
3t  ®  3z  3z 

where  T  is  the  temperature,  S  the  salinity,  the 
downward  flux  of  solar  radiation,  uw  and  vw  the  x- 
and  y-components  of  the  current  velocity,  w  the  z-com- 
ponent  of  the  current  velocity,  ug  and  the  x-  and 
y-component  of  the  geostrophic  current  velocity,  v  a 
diffusion  coefficient,  /  the  Coriolis  parameter,  qb  a 
reference  density  for  the  ocean  (1.025  g/cm3),  c  the 
specific  heat  for  the  ocean  (0.954  cal/g/°C),  t  the  time, 
and  z  the  vertical  coordinate.  Ensemble  means  are 
denoted  by  /  ),  and  primes  indicate  departure  from 

these  means.  Thus,  for  example,  the  quantity  w'S' 
represents  the  vertical  eddy  (i.e.,  turbulent)  flux  of 
salinity. 

The  terms  involving  v  in  Eqs.  (10)-(13)  represent  very 
weak  “background”  eddy  diffusion  that  exists  even 
below  the  mixed  layer.  The  value  of  v  is  set  to 
0.3  cm2/sec. 

The  vertical  eddy  fluxes  of  temperature,  salinity,  and 
momentum  are  given  by 


v'T' 


v 

_KH  37 


(14) 


h.  Boundary  conditions 

The  boundary  conditions  for  the  upper  snow/ice 
boundary  are  provided  by  the  surface  fluxes,  F,  Fn 
Fl,  and  Fr.  These  fluxes  are  used  in  the  heat  balance 
at  the  surface.  The  lower  boundary  condition  is  that 
the  bottom  of  the  ice  remains  at  the  freezing  point. 

The  upper  boundary  conditions  for  the  ocean  dej>end 
on  whether  it  is  ice-covered  or  ice-free.  When  ice- 
covered,  the  boundary  condition  for  the  temperature 
is  the  freezing  point  of  the  surface  waters.  The  silin- 
ity  boundary  condition  is  dependent  upon  whether  the 
ice  is  melting  or  freezing.  It  is  assumed  that  when  the 
ocean  freezes  all  the  salt  is  excluded  from  the  ice,  thus 
increasing  the  salinity  of  the  surface  layers.  When  the 
ice  melts,  fresh  water  is  added  to  the  surface  layer.  The 
salinity  flux  is  then  given  by 

as  ah. 

Tooo  at  (18) 

where  e,  is  the  density  of  ice  (0.92  g/cm3),  a S/1C00  is 
the  salinity  difference  per  mill,  assumed  to  be  30,  and 
A/i/A/  is  the  change  in  the  ice  thickness  with  ime. 
Thus,  an  increase  in  ice  thickness  means  the  saJniry 
flux  is  increased  at  the  surface,  and  a  decrease  n  ice 
thickness  means  a  decrease  in  the  salinity  flux  tit  the 
surface.  A  constant  fresh-water  influx  of  72  g/cm:/ 
year  is  also  assumed  (Pollard  et  al.,  1983). 

When  there  is  no  ice  the  upper  boundary  conditions 
for  the  temperature,  salinity,  and  momentum  equa¬ 
tions  are  provided  by  the  surface  fluxes.  Thus, 


[-^  *  “  g  .  fc*  * F.  * 


PoC 


(19) 


[-v's'  * v  II  Lo 


FjS, 


(20) 


4 


3u  1 

t 

V  T— V  -.  m 

Bz  Jz*0 

ax 

» 

(21) 

Po 

3v  I 

T 

v  r-v  _  = 

3z  Jz=0 

-a3L 

• 

(22) 

The  lower  boundary  conditions  for  the  ocean  are 
provided  by  holding  the  initial  temperature,  salinity, 
and  momentum  at  the  bottom  of  the  water  column 
constant. 

The  upper  boundary  condition  for  the  ice  velocity 
is  provided  by  the  wind  stress.  The  lower  boundary 
condition  is  provided  by  the  stress  of  the  water  velocity 
relative  to  the  ice  velocity.  The  stress  in  both  instances 
are  computed  using  the  drag  law  as  follows: 
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where  q0  is  the  density  of  air,  ua  and  vg  are  the  x-  and 
y-components  of  the  wind,  and  the  angle  25°  is  the 
value  assumed  for  the  turning  angle  between  the  ice 
and  the  water. 

i.  Solar  radiation  in  the  ocean 

When  the  ocean  is  ice-free,  solar  radiation  can  pene¬ 
trate  into  the  water  column  and  warm  the  surface 
layers.  The  solar  radiation  absorbed  at  each  depth  level 
is  determined  from  the  following  equation: 


aFSR(z)  *  (1  -  VFR  <e 
-0.00067z.  , 


-0.00067z, 


where  is  the  albedo  of  open  water,  and  z k  ana 
are  depth  levels  within  the  ocean. 

j.  Computing  oceanic  heat  flux 

The  oceanic  heat  flux  FB  is  computed  by  keeping 
track  of  the  change  in  the  heat  content  of  the  water 
column  and  the  heat  input  by  solar  radiation,  diffu¬ 
sion,  and  from  the  deep  ocean.  The  equation  used  to 
do  this  is  derived  by  first  integrating  Eq.  (10)  over  the 
water  column  to  get 


f=Zl  3T 
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w  here  the  term  involving  v  is  combined  with  the  w  ’  7 
term,  and  z,  =  500  m.  Evaluating  the  integral  yields 
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Assuming  that  w  at  the  surface  is  zero  and  soMng 
for  ivT|0,  the  heat  flux  at  the  surface,  Eq.  (39) 
becomes 

C1  If  *  *  ^ 

•  7 SR ^ Zj 

-  FSr'0  •  (30) 

Then,  the  oceanic  heat  flux  is 


FB  "  p°c  [W'T’I0] 


The  quantity  wT^  represents  heat  brought  into  the 
water  column  from  the  deep  ocean,  w'7",  repre¬ 
sents  the  exchange  of  heat  across  the  bottom  bound¬ 
ary  by  diffusion,  and  [FSf,\.  ~  ^SR  ■  oj  represent 
heating  by  solar  radiation. 
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3.  Case  1 — 

climatology  model  runs 

a.  Setup  of  models 

The  first  case  to  be  considered  is  a  repeat  of  one  of 
the  test  cases  from  Semtner  (1976).  In  particular,  the 
multiyear  ice  cycles,  shown  in  his  Figure  8  (Fig.  11, 
this  paper).  In  this  example  all  the  ice  melted,  which 
yielded  an  open-water  period  during  the  summer.  In 
the  fall  the  ice  reappeared  and  grew  rapidly  through 
the  winter.  The  ice  persisted  for  6  years  before  another 
open- water  period  occurred.  This  case  was  chosen  so 
that  our  simulation  of  Semtner’s  model  could  be 
verified  and  so  that  we  could  compare  the  effects  of 
the  different  treatments  of  the  mixed  layer  on  the 
growth  and  decay  of  ice.  Henceforth,  Semtner’s  model 
w  ill  be  referred  to  as  the  constant  mixed-layer  model 
(CML)  and  ours  as  the  variable  mixed-layer  model 
(VML). 

To  more  readily  compare  the  models,  many  of  the 
same  parameters  used  in  the  CML  model  are  used  in 
the  VML  model.  These  parameters  and  a  few  used  only 
in  the  VML  model  are  listed  in  Table  I.  The  values 
for  Ht  and  //,  are  the  initial  ice  thickness  and  snow 
thickness,  respectively.  The  oceanic  heat  flux,  FB ,  is 
a  constant  in  the  CML  model.  In  this  simulation,  its 
value  is  -6  kcal/cm: 'year  (-  1.929E-4  cal/cm:/sec 
or  -8.07  W.  m:).  In  the  VML  model,  FB  is  deter¬ 
mined  from  Eq.  (HJ  in  Section  2j.  The  heat  flux  from 
the  deep  ocean.  wT\,  however,  is  a  constant  and  is 
given  the  value  of  -  6  kcal/cm:  year.  The  surface 
temperature  Ts  is  initialized  to  -  2°C  in  both  models. 
The  temperature  at  the  bottom  of  the  ice  TB  is 
assumed  to  always  be  at  the  freezing  point  of  the  sur¬ 
face  water.  In  the  CML  model,  this  value  was  always 
-  2°C.  In  the  VML  model,  the  freezing  point  varies 
with  the  salinity  of  the  mixed  layer.  The  Stefan- 
Boltzmann  constant  o  is  about  2^0  higher  than  the 
standard  value.  Semtner  used  this  value  to  compare 
his  results  with  those  of  Maykut  and  Untersteiner 
(1969).  In  Section  4  the  effect  of  changing  this  value 
to  the  accepted  one  will  be  discussed.  The  value  given 
for  a ,  is  used  in  the  0-layer  ice  model  to  account  for 
some  of  the  solar  radiation  which  penetrates  into  snow- 
free  ice  (appendix,  Semtner,  1976).  The  quantities  qs 
and  qb  are  the  volumetric  heats  of  fusion  at  the  top 
and  the  bottom  of  the  ice.  Semtner  questioned  the 
validity  of  using  two  different  values  for  this  quantity, 
but  again  used  it  for  the  sake  of  comparison  with  May¬ 
kut  and  Untersteiner.  y  is  a  correction  factor  Semtner 
used,  since  the  0-layer  ice  model  does  not  keep  track 
of  solar  radiation  stored  in  the  ice  in  brine  pockets. 
Without  the  brine  pockets  the  ice  is  thinner  than  it 
should  be.  To  compensate  for  this  error,  Semtner 
increased  the  snow  and  ice  conductivities  by  the  amount 


1.065.  The  mixed-layer  depth  (MLD),  as  mentioned 
earlier,  was  30  m  in  the  CML  model.  In  the  VML 
model,  this  quantity  is  determined  by  the  dynamics  of 
the  mixed-layer.  The  parameter  SFF  represents  fresh¬ 
water  input.  The  time  step  for  the  0-layer  ice  model 
is  8  hours.  The  time  step  for  the  mixed-layer  model 
is  30  minutes. 

b.  Comparing  model  results 

Each  model  was  integrated  for  20  years.  The  ice 
thickness  and  ice  plus  snow  thickness  values  for  both 
cases  are  shown  in  Figure  12.  The  top  graph  in  this 
figure  is  the  CML  case.  The  bottom  graph  is  the  VML 
case. 

A  6-year,  no-ice  cycle,  such  as  the  one  obtained  by 
Semtner,  quickly  develops  in  the  CML  simulation.  In 
the  VML  simulation,  a  short  open  water  period  occurs 
during  year  2  and  again  in  year  4.  Then  beginning  in 
year  6,  open-water  occurs  every  year,  with  the  length 
of  the  open-water  period  increasing  each  successive 
year.  The  no-ice  cycles  for  these  cases  are  obviously 
quite  different.  These  differences  must  be  due  to  the 
differences  in  the  treatment  of  the  mixed  layer. 

In  the  CML  case,  the  only  time  the  mixed  layer  is 
used  is  when  open  water  occurs,  and  only  the 
temperature  of  the  mixed  layer  is  computed.  A  typica 
plot  of  mixed-layer  temperature  for  a  year  in  which 
open  water  occurs  is  shown  in  Figure  13.  The  watei 
temperature  increases  to  a  maximum  of  -  1.002'C 
during  the  open-water  period,  and  then  decreases  bad 
to  a  -2.0CC  value  before  ice  is  reformed.  When  the 
ocean  is  ice  covered  the  water  temperature  remains  nea 
the  freezing  point,  which  is  assumed  to  be  -2.0  C 
This  pattern  changes  little  throughout  the  simulation. 
From  Table  2,  one  can  see  that  the  maximum  mixed 
layer  temperature  varies  by  less  than  0.3:C  between 
the  open-water  periods. 

The  length  of  each  open-water  period  is  also  fairl . 
constant.  In  the  case  shown  in  Figure  13,  the  open- 
water  period  lasted  for  59.7  days.  The  average  open- 
water  period  during  the  simulation  is  54.7  days.  Eac  i 
year  in  which  open  water  occurs,  the  ice  cover  disap¬ 
pears  in  late  August  (around  Aug.  20)  and  re-forms 
in  mid-October  (around  October  15).  The  heaviest 
snowfall  in  the  simulation  occurs  during  the  months 
of  September  and  October  (Fig.  10).  Thus,  open  wan  r 
exists  during  the  time  when  the  snowfall  is  the  greates 
Snow  cover  acts  as  an  insulator  that  slows  down  the 
growth  of  ice.  With  a  very  thin  snow  cover,  which 
develops  in  the  latter  part  of  the  year,  this  insulatirg 
effect  is  greatly  reduced,  and  ice  growth  is  quite  rapid. 
From  Table  2,  it  is  seen  that  the  ice  thickness  is  the 
greatest  in  the  year  following  an  open-water  event.  The 
ice  thins  each  successive  year  until  open  water  occu's 
again.  The  snow  cover  during  these  years  is  thicker  the.* 
in  the  year  following  the  open-water  event.  The  thicker 
snow  cover  insulates  the  ice  from  the  surface  fluxes 
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Table  2.  Listing  of  ice  thicknesses,  number  of  open-water  days,  and  mixed- 
layer  temperatures  for  a  20-year  simulation  of  the  CML  model.  (Case  1) 


Year 

0  Open 

Water 

Days 

Mean  Annual 

Ice 

Thickness 

Maximum 

Ice 

Thickness 

T(1)  Max 

•c 

1 

0.0 

52.09 

94  13 

-2.000 

2 

59  7 

52.06 

89.31 

-1.002 

3 

0.0 

105.79 

159.85 

-2.015 

4 

0.0 

84  12 

124.65 

-2.016 

5 

0.0 

70.18 

110  61 

-2.015 

6 

0.0 

60.90 

101.35 

-2.015 

7 

0.0 

54.67 

95.26 

-2.015 

8 

51.0 

53.59 

91.12 

- 1 .267 

9 

0.0 

101.07 

153.75 

-2.002 

10 

0.0 

80.94 

121.44 

-2.002 

11 

0.0 

68.06 

108.49 

-2.002 

12 

0.0 

59.49 

99.96 

-2  002 

13 

0.0 

53.71 

94.33 

-2  002 

14 

54.7 

5302 

90.48 

-  1  154 

15 

0.0 

102.66 

156  33 

-2  004 

16 

00 

82  04 

122.55 

-2.004 

17 

0.0 

68.79 

109.22 

-2  004 

18 

0.0 

59.97 

100.44 

-2  004 

19 

0.0 

54.04 

94.65 

-2.004 

20 

53.3 

53.22 

90.70 

-1.197 

T(1)  Max— maximum  temperature  of  mixed  layer 


and  slows  down  the  freezing  rate;  however,  the  oceanic 
heat  flux  is  unaffected,  since  it  is  a  prescribed  constant. 
The  oceanic  heat  flux  then  gradually  reduces  the  ice 
thickness  by  melting  at  the  bottom. 

In  the  VML  simulation,  the  mixed  layer  interacts 
with  the  ice  or  the  atmosphere  every  time  step.  The 
heat,  salt,  and  momentum  fluxes  at  the  water  surface 
affect  heating  and  mixing  within  the  water  column. 
Heating  and  mixing  within  the  water  column  then 
affects  the  growth  and  decay  of  the  ice.  The  mixed- 
layer  temperature  for  year  2  from  the  VML  simula¬ 
tion  is  shown  in  Figure  14.  The  temperature  varies  from 
-  1.705°C  to  -  1.546°C  during  this  year.  The  max¬ 
imum  value  occuts  during  the  open-water  period.  The 
mixed-layer  temperature  shown  here  is  actually  the 
temperature  at  the  2.5-m  depth  level  (which  is  the  first 
level  in  the  vertical  grid),  since  the  depth  of  the  mixed- 
layer  deepens  and  shallows  in  response  to  forcing.  The 
mixed-layer  depth  throughout  year  2  is  shown  in  Fig¬ 
ure  15.  The  mixed-layer  depth  is  defined  as  the  max¬ 
imum  depth  at  which  the  salinity  has  changed  by  no 
more  than  0.1  ppt  from  its  surface  value.  Salinity  is 
used,  since  the  density  is  strongly  influenced  by  the 
salinity  in  the  Arctic  regions.  As  expected,  the  mixed 
layer,  as  seen  in  the  figure,  deepens  in  the  winter  and 
shallows  in  the  summer.  During  the  open-water  period, 
the  mixed  layer  is  only  2.5  m  deep.  Also,  open  water 
occurs  rather  late  in  the  year.  A  little  heating  occurs 


in  the  mixed  layer,  but  because  it  is  thin,  it  quickly 
cools  to  the  freezing  point  and  ice  re-forms.  The  open- 
water  period  is  short  in  this  case,  since  a  thin  layer  of 
water  is  being  heated  and  cooled.  The  open-wa  er 
period  is  only  4.7  days.  As  with  the  CML  case,  ihe 
open  water  occurs  late  in  August,  but  unlike  the  CML 
case,  ice  re-forms  quickly;  so  ice  exists  during  the  heavy 
snowfall  months.  The  resulting  snow  cover  insulates 
the  ice  and  slows  down  the  freezing  rate.  The  ice 
thickness  in  year  3  in  the  VML  case  is  quite  thin  ner 
than  in  year  3  in  the  CML  case  (maximum  of  94.04 
cm,  as  compared  with  a  maximum  of  159.89  cm).  The 
decrease  in  ice  growth  eventually  results  in  more  fre¬ 
quent  open-water  periods  later  in  the  simulation.  Snow 
cover  is  not,  however,  the  only  factor  in  the  thinner 
ice.  The  oceanic  heat  flux  is  also  important. 

The  oceanic  heat  flux  is  not  constant  in  the  VML 
simulation,  but  is  computed  from  the  change  in  the 
heat  content  of  the  water  column  as  described  in  Sec¬ 
tion  2j.  Table  3  is  a  listing  of  several  values  computed 
by  the  model  and  includes  the  amount  of  heat  input 
into  the  water  column  from  the  deep  ocean,  from  dif¬ 
fusion,  and  from  solar  heating.  The  heat  input  from 
the  deep  ocean  and  from  diffusion  vary  little  from  year 
to  year.  A  large  variability  occurs  in  the  solar  heating 
term.  This  results  in  variability  in  the  oceanic  heat  flux 
and  the  net  heating  of  the  water  column.  Solar  radia¬ 
tion  penetrates  into  the  water  column  as  described  in 
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Table  3.  Listing  of  ice  thicknesses,  number  of  open-water  days,  mixer-layer  temperatures,  and  oceanic  heating 
for  a  20-year  simulation  of  the  VML  model.  (Case  1) 


Year 

#  Open 
Water 
Days 

Mean  Annual 

Ice 

Thickness 

T(1)  Max 

•c 

T  Max 

•c 

O—O 

Diff 

cal/cm2/year 

Solar 

Fs 

1 

Net 

1 

0.0 

51  68 

-1.611 

-1.460 

5999.0 

259.0 

0.0 

-6193.4 

65 

2 

4.7 

50  08 

-1  546 

-1.381 

5998.4 

314  8 

982.2 

-  7005.5 

289  1 

3 

00 

53.31 

-1.625 

-1.449 

5997.8 

358.5 

0.0 

-6292  3 

64  0 

4 

4  3 

50  37 

-1.555 

-1  408 

5997  4 

392.9 

901.5 

-  7026  4 

265  • 

5 

0.0 

53  10 

-1.632 

-1.397 

5997.1 

419  8 

0.0 

-6338  3 

78.fi 

6 

9  3 

48.63 

- 1  572 

-1.105 

5996.8 

440.7 

2040  5 

-7884.2 

59311 

7 

24.0 

46  35 

-1.391 

-0.546 

5996  6 

456.9 

4988  6 

-  10454  0 

988  0 

8 

40  0 

43  66 

-0.976 

-0.123 

5996.5 

469  3 

7950.7 

-13369  0 

1047  ij 

9 

35.7 

45  26 

-1.090 

-  0  224 

5996  4 

478  5 

7176  5 

-13103  0 

548  i 

10 

38  3 

44  68 

- 1  020 

-0.156 

5996  3 

480  5 

7648  9 

-13526  0 

604  ) 

11 

41  3 

43  84 

-0  924 

-0.065 

5996.3 

489.3 

8218  7 

-14082  0 

622  > 

12 

43  0 

43.11 

-0.850 

0  005 

5996  2 

491.5 

8591.3 

-  14482  0 

597  3 

13 

46.7 

42.11 

-0  724 

0  123 

5996  2 

492,0 

9294  7 

-15145  0 

638  3 

14 

49  0 

41.36 

-0.619 

0  221 

5996  2 

491.0 

9807.7 

-  15659  0 

636  3 

15 

52  3 

39  99 

-0.463 

0.363 

5996  3 

488.9 

10573.0 

-16419  0 

639 .3 

16 

54  3 

39  00 

-0.353 

0  464 

5996  3 

485  8 

11073.0 

-16934  0 

621  5 

17 

53  0 

37  72 

-0.298 

0  517 

5996.3 

481.8 

11173.0 

-16838  0 

813  5 

18 

62  0 

33  15 

0.161 

0  946 

5996  4 

477.2 

13347.0 

-18552  0 

1269  D 

19 

67  3 

30  56 

0  458 

1.227 

5996  5 

472  0 

14732.0 

-  19978  0 

1222  9 

20 

80  7 

30  04 

0  672 

1.431 

5996  5 

466  3 

16077.0 

-  22374  0 

165  8 

T(1)  Max— maximum  temperature  at  depth  level  1  (2  5  m) 

T  Max— maximum  temperature  in  the  water  column 

How,c — '°tal  heat  input  into  the  water  column  py  advection  from  the  deep  ocean 
Diff— total  heat  input  into  the  water  column  by  diffusion  across  the  bottom 
Solar— total  heat  input  into  the  water  column  by  solar  radiation 
FB— total  heat  transfer  across  the  ice/ocean  or  ice/water  interface 
Net— +  Diff  +  Solar  +  Fb 


Section  2i.  Each  depth  level  absorbs  different  amounts 
of  heat  from  solar  radiation.  The  surface  level  (2.5-m 
level)  absorbs  the  most  solar  heat.  It  also  reacts  the 
fastest  with  the  atmospheric  forcing;  thus,  the  heating 
or  cooling  of  the  surface  layer  can  be  quite  rapid. 
Figure  16  shows  some  temperature  profiles  for  steps 
before  (step  695),  during  (steps  700-720),  and  after 
(steps  725  and  730)  the  open-water  period.  One  sees 
that  the  water  column  is  being  warmed  during  the  open- 
water  period.  The  subsurface  waters  are  warmed  the 
most,  since  they  are  slower  to  respond  to  surface 
forcing.  The  subsurface  layers  are  insulated  from  the 
surface  forcing  at  this  time  because  the  water  column 
is  highly  stratified  and  because  no  mixing  occurs  until 
ice  starts  to  grow  again  (Fig  17).  The  maximum 
temperature  that  occurs  in  the  water  column  in  this 
case  is  -  1.381°C  and  occurs  at  the  12.5-m  depth  level 
at  time  step  720. 

Ice  begins  to  grown  again  at  step  721,  since  the 
temperature  in  the  surface  waters  has  dropped  below 
the  freezing  point.  The  temperature  in  the  subsurface 
waters,  however,  remains  slightly  above  the  freezing 


point  after  ice  has  re-formed  (steps  725  and  730).  "This 
heat  in  the  subsurface  waters  is  released  as  the  mixed 
layer  deepens  as  the  ice  grows.  This  can  be  seen  in 
Figure  18,  which  shows  some  temperature  profiles  for 
the  last  part  of  year  2.  During  this  time  ice  grows  f  om 
1  cm  to  approximately  40  cm  thick  and  the  mixed  layer 
deepens  from  2.5  m  to  32.0  m.  The  release  of  this  heat 
in  the  subsurface  waters  during  the  growth  of  ice  slows 
down  the  freezing  rate  of  ice.  This  effect  yields  thin¬ 
ner  ice,  and  in  combination  with  the  snow  cover, 
explains  why  years  2  and  3  are  so  different  between 
the  CML  and  the  VML  simulation. 

Year  4  in  the  VivlL  simulation  is  quite  similar  to 
year  2.  In  year  6,  however,  open  water  occurs  slightly 
earlier  in  the  year  than  in  year  2  or  4,  and  it  lasts  a 
few  days  longer.  The  warming  of  the  water  colimn 
is  thus  greater,  and  the  ice  thickness  for  this  year  is 
thinner  than  in  the  previous  years.  Following  this  year 
open  water  occurs  every  year.  The  length  of  the  open- 
water  period  tends  to  increase  each  successive  year,  and 
so  does  the  heating  in  the  water  column. 
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Using  year  10,  as  an  example,  the  heating  of  the 
water  column  will  be  examined  in  more  detail.  The 
open-water  period  this  year  lasted  38.3  days.  It  opened 
up  on  August  10  and  began  to  refreeze  on  Septem¬ 
ber  20.  This  period  is  not  38.3  consecutive  days  of  open 
water,  but  36  days  of  open  water  followed  by  a  brief 
period  of  ice  growth  between  September  16  and 
September  18.  During  this  time  2  cm  of  ice  formed. 
The  mixed  layer  deepened  as  a  result  of  the  ice  growth 
releasing  some  heat  from  below  the  mixed  layer.  The 
heat  was  sufficient  to  melt  the  thin  ice  cover  and  to 
keep  the  ocean  ice-free  until  September  20  when  con¬ 
tinuous  ice  growth  began. 

Temperature  profile  plots  for  this  year  are  shown 
in  Figures  19  and  20.  At  step  655  in  Figure  19,  no  warm 
spot  is  observed  in  the  water  column  down  to  100  m. 
At  step  680,  the  warm  spot  begins  to  form  and  grows 
until  step  780  (open  water  is  from  step  659  to  step  780). 
After  this  step,  the  warm  spot  is  seen  to  decrease  at 
the  same  time  that  the  mixed  layer  deepens.  However, 
not  all  the  heat  in  the  warm  spot  is  exhausted  by  the 
end  of  the  year  (step  1080  in  Fig.  20).  In  fact,  it  is  not 
until  early  spring  of  the  following  year  that  the  warm 
spot  disappears  (Fig.  21:  step  375  occurs  on  May  6). 
Thus,  the  heat  input  into  the  water  column  during  this 
open-water  period  affects  the  growth  of  ice  for  several 
months  after  ice  re-forms.  The  net  result  is  thinner  ice 
and  warmer  water  below  the  ice.  In  general,  as  seen 
in  Table  3,  the  heat  input  into  the  water  column  in  a 
year’s  time  is  not  all  released  in  that  year.  The  previous 
year  greatly  affects  the  following  year. 

Examining  Table  3,  one  sees  that  the  solar  heat  input 
tends  to  increase  with  each  successive  year,  beginning 
with  year  9,  and  that  the  oceanic  heat  flux,  FB,  also 
tends  to  increase,  except  between  years  16  and  17.  The 
net  heating,  however,  fluctuates  from  year  to  year.  This 
implies  that  the  release  of  heat  from  the  water  column 
is  not  in  a  steady  manner,  but  is  highly  variable. 

c.  Oceanic  heat  flux 

The  variability  of  the  oceanic  heat  flux  throughout 
a  year  is  shown  in  Figure  22.  The  dashed  line  in  the 
plot  is  the  value  -6.0  kcal/cm2/year,  the  constant 
heat  input  from  the  deep  ocean.  The  oceanic  heat  flux 
throughout  much  of  this  year  is  considerably  larger 
than  this  constant  value.  The  only  time  FB  remains 
near  this  constant  value  is  in  the  early  summer  when 
ice  is  melting  and  the  mixed-layer  is  shallowing.  The 
largest  values  of  oceanic  heat  flux  occur  in  the  winter 
following  the  open-water  period.  Comparing  this  plot 
to  the  plot  of  the  salinity  mixed  layer  depth  (Fig.  23), 
one  sees  that  these  large  oceanic  heat  flux  values  cor¬ 
respond  to  times  when  the  mixed  layer  is  deepening. 
So,  just  as  was  seen  in  the  temperature  profiles,  the 
heat  stored  in  the  water  column  is  released  as  the  mixed 
layer  deepens.  From  Figure  22,  it  is  clear  that  this  heat 
is  released  in  bursts  and  that  some  of  these  bursts  are 
quite  large. 


The  features  shown  in  year  10  are  repeated  in  other 
years.  Figure  24  is  a  plot  of  the  last  5  years  of  the 
simulation  showing  the  oceanic  heat  flux,  the  salinity 
mixed-layer  depth,  and  the  ice  thickness.  The  main  d  f- 
ferences  between  these  years  is  in  the  length  of  the  open- 
water  periods,  the  deepening  of  the  mixed  layer,  and 
the  magnitude  of  the  oceanic  heat  flux  bursts.  Oie 
interesting  feature  is  the  spikes  in  the  mixed  layer  depth 
plot.  These  spikes  correspond  to  times  in  which  tie 
mixing  process  brings  up  sufficient  heat  to  melt  a  few 
centimeters  of  ice.  A  thin  layer  of  low-salinity  water 
thus  forms  in  the  surface  waters.  The  mixed-layer  depth 
then  shallows  briefly.  The  continued  growth  of  ice  so  jn 
results  in  mixing  of  the  surface  waters,  which  brings 
the  mixed-layer  depth  back  to  its  previous  depth. 

In  the  VML  model,  the  deepening  of  the  mixed  layer 
is  controlled  by  density  instabilities  and  velocity  shears 
in  the  water  column.  The  density  instabilities  are  cc  n- 
trolled  primarily  by  the  salinity  flux,  since  in  the  Arctic 
the  temperature  usually  changes  less  than  the  salini  y. 
In  this  simulation,  when  ice  is  growing  the  salinity  and 
thus  density  of  the  surface  layers  increase.  Eventually 
the  density  structure  becomes  unstable  and  mixing 
occurs.  If  warmer  water  exists  below  the  mixed  layer, 
this  water  is  mixed  into  the  mixed  layer,  and  the  heat 
in  the  water  tries  to  increase  the  temperature  of  i  he 
mixed  layer.  The  temperature  of  the  mixed  layer  is  con¬ 
strained  when  ice  cover  exists,  since  the  temperature 
at  the  bottom  of  the  ice  must  remain  at  the  freez  ng 
point.  So  this  heat  either  slows  down  the  freezing  rate 
or  melts  a  few  centimeters  of  ice.  The  mixed  layer  then 
stabilizes  until  density  instabilities  arise  again  from  :he 
salinity  flux. 

Velocity  shears,  which  arise  from  the  relative  mot  on 
of  the  ice  and  water,  also  generate  mixing  in  the  water 
column.  In  this  simulation,  the  velocity  shears  are  tot 
strong  enough  to  overcome  the  density  structure,  aut 
can  greatly  affect  the  rate  of  mixing  and  the  ratt  at 
which  the  heat  is  released  from  the  water  column.  This 
is  evident  from  Figure  25,  which  is  a  plot  of  the  ocea  nic 
heat  flux  for  year  10  for  a  simulation  in  which  the 
winds  were  shut  off;  no  winds,  thus  no  velocity  shears 
occur  in  the  water  column.  This  plot  is  similar  to  Fig¬ 
ure  22,  except  for  the  spikes  observed  in  the  latter.  The 
oceanic  heat  flux  values  in  Figure  25  appear  to  vary 
more  gradually  and  to  have  fewer  jumps. 

Table  4  is  a  listing  of  some  important  quantities  fiom 
the  no-winds  simulation.  Things  change  when  the 
velocity  shears  are  removed  from  the  simulation.  The 
most  noticeable  changes  are  that  the  net  heating  va.ues 
are  generally  larger,  the  open-water  periods  tend  to  be 
longer,  and  the  water  temperature  tends  to  be  warmer 
in  the  no-winds  case.  Starting  at  year  12  in  the  no-w  nds 
case,  a  3-year  cycle  is  observed  in  the  net  heating,  the 
number  of  open-water  days,  and  the  temperature  of 
the  water  column.  The  net  heating  decreases  from  a 
large  positive  value  to  a  negative  value  over  the  3-year 
cycle.  The  number  of  open-water  days  increases  and 
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Table  4.  Listing  of  ice  thicknesses,  number  of  open-water  days,  mixed-layer  temperatures,  and  oceanic  heating 
for  a  20-year  simulation  of  the  VML  model  when  no  shears  are  developed  in  the  water  column  (Case  1 D). 


#  Open 
Water 

Mean  Annual 
Ice 

T(1)  Max 

T  Max 

cal/cm*/year 

Year 

Days 

Thickness 

°C 

°C 

Diff 

Solar 

PB 

Net 

1  0.0  52.80  -1.612  -1.460  5999.0  259.8  0.0  -  6097.5  161.3 

2  3  7  51.06  -1.563  -1.437  5998.4  314.8  754.3  -  6776.1  291.3 


3 

0  0 

55.32 

-1.628 

-1  446 

5997.8 

3585 

0.0 

-6292.5 

63.7 

4 

00 

50.57 

-1.628 

-1.419 

5997.4 

392.8 

0.0 

-6213.4 

176.8 

5 

17.7 

47.77 

-1.383 

-0.759 

5997.1 

419.7 

3719.3 

-8999.1 

1137,0 

6 

30.0 

46.27 

-1.142 

-0.351 

5996.8 

440.6 

6185.1 

-11394.0 

1228,2 

7 

40.3 

43.44 

-0.843 

0012 

5996.6 

456.6 

8355.4 

-13444.0 

1365  0 

8 

49  0 

40.54 

-  0.474 

0.360 

5996.5 

468.7 

10242.0 

-15264.0 

1443  0 

9 

55  3 

38  07 

-0.165 

0.652 

5996  4 

477.4 

11715.0 

-16744.0 

1445  0 

10 

61  0 

35  77 

0,127 

0  930 

5996.3 

482.9 

13079.0 

-  18082.0 

1476  3 

11 

65  7 

33  84 

0  366 

1.161 

5996.3 

485.7 

14201.0 

-19241.0 

1442  1 

12 

69  0 

32  37 

0.553 

1.340 

5996.3 

486.0 

15046.0 

-20161.0 

1368  0 

13 

81  3 

32  61 

0.710 

1,492 

5996.3 

484.1 

16092.0 

-22139.0 

433  3 

14 

54  3 

40  55 

-0  211 

0.626 

5996  4 

480.1 

11456.0 

-18097.0 

-164.8 

15 

66  3 

33  72 

0.402 

1.200 

5996  4 

474.2 

14347.0 

-19311  0 

1506.7 

16 

81  3 

32  79 

0  709 

1.493 

5996  5 

467.2 

16092.0 

-22080.0 

475  3 

17 

55  7 

40  19 

-0  154 

0.680 

59966 

459  1 

11739.0 

-18365  0 

-1703 

18 

67  3 

33  26 

0  464 

1.260 

5996  7 

450.4 

14622.0 

-19607  0 

1462  4 

19 

83.0 

32  15 

0.773 

1.556 

5996.8 

441.3 

16379  0 

-22515  0 

302  0 

20 

55  0 

40  46 

-0.181 

0  656 

5996.9 

432  1 

11598  0 

-18365.0 

-338  8 

T(1)  Max— maximum  temperature  at  depth  level  1  (2  5  mj 
T  Max— maximum  temperature  in  the  water  column 

H  —  total  heat  input  into  the  water  column  by  advection  from  the  deep  ocean 
Diff— total  heat  input  into  the  water  column  by  diffusion  across  the  bottom 
Solar— total  heat  input  into  the  water  column  by  solar  radiation 
Fb— total  heat  transfer  across  the  ice/ocean  or  ice/water  interlace 
Net— +  Ditl  +  Solar  +  Fe 

then  decreases  during  this  time.  It  varies  from  about  greater  than  the  change  in  temperature  in  the  water 

55  to  81  days.  The  water  temperature  also  increases  column.  In  general,  these  profiles  show  that  sound 

and  then  decreases.  This  value  varies  from  about  energy  is  trapped  in  a  duct  near  the  surface.  A  sound 
-  0.2°C  to  about  0.7°C.  channel,  however,  forms  around  50  m  depth  when  the 

water  column  has  been  warmed  during  an  open-water 
d.  Sound  Speed  event.  (Wilson’s  (I960)  equation  was  used  to  compute 

The  discussion  of  the  results  for  Case  1  will  be  con-  sound  speed.) 
eluded  with  a  brief  examination  of  sound  speed  pro¬ 
file  plots  for  year  20  of  the  simulation.  These  plots  are  A  „  ,  .  _ 

shown  in  Figure  26.  The  corresponding  temperature  4.  til  ICC  t  OI  CtlHIlglllg  O  3110  Qj) 
and  salinity  profiles  are  shown  in  Figures  27  and  28.  It  was  mentioned  earlier  that  the  values  of  o  and  qh 

These  profiles  were  chosen  to  depict  the  differing  were  not  the  generally  accepted  values.  They  were  used 

physical  conditions  throughout  the  year.  The  profiles  by  Semtner  for  the  sake  of  comparison  with  Maykut 

at  steps  1  thru  300  and  steps  900  and  1050  were  taken  and  Untersteiner.  The  value  used  for  o  is  2 °7o  higher 

during  ice  growth.  The  profile  at  step  450  was  taken  than  the  standard  value.  This  increased  value  causes 

during  ice  melt,  and  those  at  steps  600  and  750  were  more  heat  in  the  form  of  longwave  radiation  to  be 

taken  during  open  water  (open  water  lasts  from  step  released  from  the  snow/ice  system,  resulting  in  thicker 

598  to  840).  ice  than  if  the  standard  value  had  been  used.  The 

The  sound  speed  in  these  profiles  is  clearly  domi-  value  used  for  qb  was  questioned  by  Semtner.  Fie 

nated  by  the  temperature.  Some  salinity  influence,  stated  that  the  reduced  value  of  qb  at  the  bottom  of 

however,  is  observed  at  the  surface  in  the  profile  at  the  ice  “implied  an  additional  loss  of  heat  over  what 

step  600.  Here  the  decrease  in  salinity  due  to  melting  is  is  prescribed  in  surface  fluxes  must  occur  before  the 
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ice  reaches  the  upper  surface.”  In  later  work,  he  used 
the  value  of  72  cal/craJ  for  both  the  upper  and  lower 
boundaries  of  the  ice,  and  the  standard  value  of  a  (per¬ 
sonal  communication).  In  this  section,  the  effects  on 
the  VML  model  of  changing  these  parameters  will  be 
examined.  Table  5  is  a  listing  of  the  mean  annual  ice 
thicknesses  and  the  number  of  open-water  days  for  four 
cases:  CASE  1 ,  the  standard  case  as  described  earlier; 
CASE  1A,  where  qb  is  the  same  as  qs;  CASE  IB, 
where  a  is  changed  to  the  standard  value  of  1 .3545E-12 
cal/cm2/sec/0K4  and  CASE  1C,  where  both  o  and  qb 
are  changed. 

Both  reducing  o  by  2%  and  changing  qb  to  72 
causes  more  heat  to  remain  in  the  snow/ice  system, 
which  results  in  thinner  ice.  In  this  example,  it  results 
in  an  average  decrease  in  the  mean  annual  ice  thickness 
of  17.3%  for  Case  1  A,  23.7%  for  Case  IB,  and  26.6% 
for  Case  1C.  The  no-ice  cycle  also  changes  signifi¬ 
cantly.  Each  subsequent  case  has  more  frequent  and 
longer  open  water  periods.  (In  the  rest  of  the  report, 
any  model  runs  discussed  will  have  the  corrected  values 
for  a  and  qb.) 

5.  Deep  ocean  heat  flux  (Hdeep) 

Open  water  that  occurs  every  year  is  reasonable  for 
the  marginal  ice  zones,  since  the  deep  ocean  heat  flux, 
Hdnp,  can  be  quite  large.  The  exact  value  of  this  heat 
flux  depends  on  what  Arctic  area  is  being  considered. 
It  can  vary  from  near  zero  in  the  central  Arctic 
(Aagaard,  1981;  McPhee  and  Untersteiner,  1982)  to 
several  hundred  W/m2  in  the  Greenland  Sea  (Hibler 
and  Bryan,  1984).  At  this  point  in  the  study,  we  do 
not  wish  to  simulate  a  specific  region,  but  to  examine 
how  changes  in  the  prescribed  deep  ocean  heat  flux 
can  affect  the  ice/ocean  interaction. 

In  Figure  29  is  an  ice  and  ice  plus  snow  thickness 
plot,  and  an  oceanic  heat  flux  plot  for  a  20-year  simula¬ 
tion  of  the  VML  model  using  1.5  kcal/cm2/year  (2.02 
W/m:)  as  the  deep  ocean  heat  input  value.  One  can 
see  that  heat  in  the  water  column  is  not  sufficient  in 
this  case  to  cause  open  water.  In  fact,  the  mean  annual 
ice  thickness  increases  until  it  reaches  a  value  of 
115.8  cm  in  year  20.  Spikes  occur  in  the  oceanic  heat 
flux-plot  even  though  there  are  no  open-water  periods. 
This  is  because  the  mixed  layer  always  deepens  in  the 
winter,  and  in  this  simulation  the  temperature,  in 
general,  increases  with  depth.  So  warmer  water  is 
always  below  the  mixed  layer.  Figures  30  through  32 
are  similar  plots  for  simulations  using  3.0  kcal/cm2/ 
year  (4.04  W/m2),  4.0  kcal/cm2/year  (5.38  W/m2), 
and  6.0  kcal/cm2/year  (8.07  W/m2).  It  is  obvious 
that  as  the  value  of  Hdnp  increases,  the  frequency  of 
the  open- water  periods  also  increases.  Though  not  as 
obvious  the  mean  annual  ice  thickness  decreases  as 
Hdetp  increase,  and  the  length  of  the  open-water 
period  increases.  In  all  the  runs,  the  oceanic  heat  flux 


FB  hovers  near  the  deep  ocean  heat  flux  value,  Hdetp, 
except  for  the  burst  of  heat  that  occurs  when  the  mixed 
layer  deepens  and  vents  heat  stored  in  the  water 
column.  These  bursts  become  more  frequent  and  larger 
as  the  value  of  H^r  increases,  and  open  water 
becomes  more  frequent. 

6.  Case  2 — solar  radiation 
stored  in  brine  pockets 
a.  Setup  of  model 

The  variable  mixed-layer  model  is  extended  by  simu- 
lating  the  heat  stored  in  the  ice  in  brine  pockets  when 
solar  radiation  penetrates  into  snow-free  ice.  To  do 
this,  a  heat  reservoir  was  created.  Heat  is  added  to  this 
reservoir  whenever  the  ice  is  snow-free,  and  the  sur¬ 
face  temperature  of  the  ice  is  above  freezing.  The 
following  equation  is  used  to  determine  the  amount 
of  heat  stored  in  the  reservoir  each  time  step  w  hen  :he 
ice  is  snow-free: 

fBR-  “'V  '•  F»  «  -  e'°-I5h‘>'»  .I32| 

where  I0  is  the  percentage  (17%)  of  solar  radiation 
that  penetrates  into  the  ice.  It  is  also  assumed  that  some 
of  the  penetrating  solar  radiation  goes  through  the  ice 
and  enters  the  water  column.  The  amount  of  heat  enter¬ 
ing  the  water  column  is  determined  from 

FR„-  <1  -  I,  Fr  e-°-15  hi  it  .  ,33) 

Thus  for  thick  ice  most  of  the  penetrating  solar  radia¬ 
tion  is  stored  in  the  brine  pockets,  but  very  little  radia¬ 
tion,  if  any,  enters  the  water  column.  For  thin  ice, 
however,  solar  radiation  warms  the  water  column  and 
is  stored  in  brine  pockets.  This  heat  reservoir  remits 
in  less  heat  being  available  for  melting  at  the  surface 
during  the  summer.  It  can  also  delay  the  growth  of 
ice  in  the  fall  until  all  the  heat  in  the  reserve  ir  is 
released. 

Some  limitations  were  put  on  the  heat  reservoir. 
First,  heat  is  not  allowed  to  accumulate  in  the  reser¬ 
voir  when  the  ice  thickness  is  less  than  25  cm.  Second, 
a  maximum  value  is  set  for  the  amount  of  heat  that 
can  be  stored  in  the  brine  pockets.  This  maximum 
depends  on  the  thickness  of  the  ice  and  is  determined 
as  follows: 

H  =  0.5  (h.  -  h  )o 
max  i  l  min^  qb  ’  (34) 

where  himjn  is  the  minimum  ice  thickness  (25  cm)  for 
which  solar  radiation  accumulates  in  the  reservoir. 
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Table  5.  Listing  of  ice  thicknesses  and  number  of 
open  water  days  for  a  10-year  simulation  of  the  VML 
model  when  the  values  of  a  and  qb  are  changed. 


Year 

Mean  Annual  tee  Thtckneea  (cm) 
Case  1  Case  1 A  Case  IB 

Case  1C 

1 

51.68 

49  63 

43.67 

42.04 

2 

50.08 

44.06 

36.92 

34.01 

3 

53  31 

41 .80 

40.32 

33.90 

4 

50.37 

39  26 

40.77 

39.14 

S 

53.10 

40.20 

31.59 

29.51 

6 

48  63 

39.57 

39.44 

36  78 

7 

46  35 

38.61 

35.72 

34.61 

6 

43.66 

37.61 

33.06 

38.31 

9 

45  26 

36.76 

37.28 

34.16 

10 

44  68 

35  43 

32.87 

34.96 

Average 

48  71 

40.29 

37.16 

35.74 

Number  of  Open-Water  Oays 

Year 

Case  1 

Case  1A 

Case  18 

Case  1C 

1 

00 

00 

56  3 

57.7 

2 

4.7 

190 

76  3 

780 

3 

00 

32  7 

70,7 

890 

4 

4  3 

46  3 

67  3 

680 

5 

00 

45  3 

95  3 

96  7 

6 

9  3 

47  0 

73  3 

76  0 

7 

24  0 

49  3 

800 

88. 3 

8 

40  0 

523 

930 

73  7 

9 

35  7 

54  3 

77  7 

880 

10 

383 

53  7 

93  7 

860 

Table  6.  Listing  of  parameters  and  assigned  values 
used  in  Case  2  for  the  VML  model.  Brine  pockets 
and  penetrating  solar  radiation  are  included  in  this 
simulation. 


Variable  Mixed  Layer 

H 

50  cm 

H 

t 

30  cm 

Fb 

-  9.6451  E-5  ±  heal  change  (cal/cm2/sec) 

Ts 

-  2.0°C 

Ts 

Varies  with  salinity 

0 

1  3545E-12  cal/cm*/°K4/sec 

a 

0.64 

a 

w 

0.07 

q. 

72  cal/cm3 

72  cal/cm3 

r 

1.0 

MLD 

Variable 

SFF 

3.341  E-5  cm/sec 

Time  Step 

8  hours/30  minutes 

Forcing 

Heat  Fluxes: 

Climatology 

Winds 

FNOC  average  monthly 

T&S  Profiles 

AIDJEX  data  (Beaufort  Sea) 

b.  Results 


Case  1  —o  *  1  3850E-12  and  q0  -  64.0 

Case  1A— o  =  1  3850E-12  and  qo  -  72  0 

Case  iB—o  *  1  3545E-12  and  q6  -  64.0 

Case  1C— o  -  1  3545E- 1 2  and  qo  =  72  0 

The  heat  in  the  reservoir  is  released  whenever  the 
surface  temperature  of  the  ice  falls  below  freezing.  Just 
enough  heat  is  released  to  raise  the  surface  temperature 
to  the  freezing  point.  Thus  the  surface  temperature  of 
the  ice  is  maintained  at  the  freezing  point  until  all  the 
heat  in  the  reservoir  is  released.  For  thin  ice  the  max¬ 
imum  allowed  heat  storage  is  usually  reached  before 
the  surface  temperature  falls  below  freezing.  When  this 
occurs,  no  me  -  heat  is  diverted  to  the  heat  reservoir 
(/0  =  0.0),  and  heat  in  the  reservoir  is  gradually 
released.  This  heat  is  treated  as  additional  ice  melt.  This 
simulation  is  similar  to  that  used  by  Semtner  in  his 
model. 

Two  of  the  parameters  listed  in  Table  1  had  to  be 
changed  for  the  simulation  of  the  storage  of  solar  radia¬ 
tion  in  brine  pockets.  These  parameters  are  ».  and  y. 
The  value  of  tr,  is  changed  to  0.64  and  y  is  changed 
to  1 .0.  Both  of  these  parameters  had  larger  values  in 
Case  1  to  loosely  simulate  solar  radiation  being  stored 
in  brine  pockets.  (See  Table  6  for  a  list  of  all  the 
parameters  used  in  this  run.) 


Thi$  version  of  the  VML  model  (Case  2)  was  inte¬ 
grated  for  20  years.  The  results  of  this  simulation  will 
be  compared  with  those  for  the  case  shown  in  Figure  30 
(Case  IF).  The  ice  and  ice  plus  snow  thickness  plot  for 
Case  2  is  shown  in  Figure  33.  Comparing  Figures  30 
and  33,  it  is  apparent  that  thinner  ice  and  more  fre¬ 
quent  open-water  periods  occur  when  penetrating  solar 
radiation  is  simulated  than  when  it  is  not.  Tables  7  and 
8  list  several  important  values  obtained  from  the 
simulations.  In  Case  2  open-water  occurs  beginning 
at  year  1.  By  year  11,  the  length  of  the  open-water 
period  each  year  alternates  between  approximately  88 
and  67  days.  In  case  IF,  open  water  does  not  occur 
until  year  20,  and  then  for  only  two  days.  The  mean 
annual  ice  thicknesses  are  on  the  order  of  10  to  20  cm 
thinner  in  Case  2  than  in  Case  IF,  and  the  water 
temperatures  are  a  few  degrees  warmer  in  Case  2  than 
in  Case  IF.  The  reason  for  these  differences  is  seen 
in  the  heat  flux  values  listed  in  the  tables.  The  con¬ 
trolling  difference  is  the  solar  values.  In  Case  IF,  no 
solar  heating  of  the  water  column  occurs  until  year  20. 
In  Case  2,  however,  solar  heating  is  quite  pronounced 
and  begins  at  year  l .  Several  thousands  of  cal/cm2  of 
heat  are  input  into  the  water  column  in  a  year’s  time 
by  solar  heating.  This  is  particularly  countered  by  an 
increase  in  the  oceanic  heat  flux,  FB.  The  net  heating 
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Table  7.  Listing  of  ice  thicknesses,  number  of  open-water  days,  mixed-layer  temperatures,  and  oceanic  heating 
for  a  20-year  simulation  of  the  VML  model  when  no  brine  pockets  or  penetrating  solar  radiation  is  allowed 
(Case  IF). 


Year 

#  Open 
Water 
Days 

Mean  Annual 
Ice 

Thickness 

T(1)  Max 
°C 

T  Max 

«c 

Ditt 

cal/cm2/year 

Solar 

Fe 

Net 

1 

0.0 

55.16 

-1.611 

-1.460 

2999.5 

260.1 

0.0 

-3226.1 

33.5 

2 

0.0 

57.15 

-1.619 

-1.487 

2999.2 

315.7 

0.0 

-3195.7 

1  "9.2 

3 

0.0 

58.32 

- 1 .623 

-1.458 

2996.9 

360.0 

0.0 

-3225.0 

133.9 

4 

0.0 

58.80 

- 1 .626 

-1.432 

2998.7 

395.0 

0.0 

-3254.7 

139.0 

5 

0.0 

58.78 

-1.628 

-1.410 

2998.5 

422.4 

0.0 

-3284.0 

136.9 

6 

0.0 

58.41 

-1.629 

-1.391 

2998.4 

443.8 

0.0 

-3311.7 

130.5 

7 

0.0 

57.84 

-1  630 

- 1 .376 

2998.3 

460.5 

0.0 

-3336.2 

122.5 

8 

0.0 

57.17 

-1.631 

-1.363 

2998.2 

473.4 

O.C 

-3357  2 

114.4 

9 

0.0 

56.47 

-1.631 

- 1 .352 

2998.2 

483.3 

0.0 

-3375.0 

106.4 

10 

0.0 

55.78 

-1.632 

-1.343 

2996.1 

490.8 

0.0 

-3389  9 

99.0 

11 

0.0 

55.12 

-1.632 

-  1.336 

2998.1 

496.5 

0.0 

-3402  4 

92.1 

12 

0.0 

54.52 

-1.632 

-  1.329 

2998.1 

500.6 

0.0 

-3412.9 

85.8 

13 

0.0 

53.98 

-1.632 

-1.324 

2998.0 

503.6 

0.0 

-3421.6 

80.1 

14 

0.0 

53  49 

-1.632 

-1.319 

2998.0 

505  7 

0.0 

-3428  9 

74.9 

15 

0.0 

53  06 

-1.632 

-1.315 

2998  0 

507.1 

0.0 

-3435.0 

70.1 

16 

0.0 

52.69 

-1.632 

-1.312 

2998.0 

507.9 

0.0 

-3440.1 

65.8 

17 

0.0 

52.36 

-1.631 

-  1.309 

2998.0 

508.3 

0  0 

-3444.5 

61.9 

18 

00 

52.07 

-1.631 

-  1.306 

2998.0 

508.4 

0.0 

-3448  1 

58.3 

19 

00 

51.83 

- 1.631 

-1.304 

2998  0 

508.2 

0.0 

-3451.3 

550 

20 

20 

52  85 

-1.630 

-1.300 

2998  0 

507  9 

390.7 

-3851  5 

45.1 

T(1)  Max— maximum  temperature  at  depth  level  1  (2.5  m) 

T  Max— maximum  temperature  in  the  water  column 

— total  heat  input  into  the  water  column  by  advection  from  the  deep  ocean 
Ditt— total  heat  input  into  the  water  column  by  diffusion  across  the  bottom 
Solar— total  heat  input  into  the  water  column  by  solar  radiation 
FB— total  heat  transfer  across  the  ice/ocean  or  ice/water  interface 
Net— +  Diff  +  Solar  +  FB 


of  the  water  column  and  the  temperature  of  the  water  water  column  is  apparent  from  Figure  34,  which  plots 

column  are  thus  greater  in  Case  2  than  in  Case  IF.  temperature  profiles  for  Case  IF  (the  dashed  line)  and 

Year  1  of  the  above  simulation  will  be  examined  in  Case  2  (the  solid  line)  at  several  time  steps  during  the 

greater  detail  to  get  a  better  understanding  of  how  the  ice-melting  phase.  At  step  502  the  profiles  are  iden- 

penetrating  solar  radiation  simulation  affects  the  ice  tical.  At  later  time  steps,  the  temperature  of  the  sur- 

and  the  ocean.  The  simulations  of  Case  IF  and  Case  2  face  waters  for  Case  2  is  seen  to  increase  due  to  solar 

are  nearly  identical  for  the  first  part  of  year  1 .  A  slight  heating.  The  increase  is  not  great,  but  it  does  have  an 

difference  in  ice  thicknesses  results  due  to  the  difference  effect  on  the  growth  of  ice. 

in  the  value  of  y.  This  parameter  affects  the  conduc-  The  heat  stored  in  the  heat  reservoir  gradually 

tion  of  heat  in  the  ice.  increases  until  a  maximum  value  is  reached.  A  max- 

The  snow  is  totally  melted  by  step  502  in  both  cases.  imum  value  tied  to  the  thickness  of  the  ice  is  set  for 

At  this  point  the  penetrating  solar  radiation  Simula-  the  heat  reservoir,  so  that  the  heat  in  the  reservoir  can- 

tion  begins  to  act.  Table  9  lists  the  amount  of  ice  melt  not  exceed  the  amount  needed  to  melt  all  the  ice.  In 

and  the  ice  thickness  at  several  time  steps  during  the  this  case  a  maximum  value  of  644.6  cal/cm-  for  the 

melting  process.  The  amount  of  ice  melted  at  the  sur-  heat  reservoir  is  reached  at  step  595.  After  this  step, 

face  is  decreased,  and  that  melted  at  the  bottom  is  the  heat  in  the  reservoir  decreases  and  the  total  ice  melt 

increased  in  Case  2  as  was  expected.  The  total  melt,  increases.  The  ice  melt  listed  in  Table  9  from  step  602 

however,  is  less  in  Case  2  than  in  Case  IF  because  some  on  shows  some  increase  in  the  melting  at  the  bottom 

heat  is  stored  in  the  heat  reservoir,  and  some  goes  to  and  at  the  surface  of  the  ice  for  Case  2.  The  values 

warming  of  the  water  column.  The  warming  of  the  listed  in  the  table  are  only  the  melt  due  to  the  heat 
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Table  8.  Listing  of  ice  thicknesses,  number  of  open-water  days,  mixed-layer  temperatures,  and  oceanic  heating 
for  a  20-year  simulation  of  the  VML  model  when  brine  pockets  and  penetrating  solar  radiation  are  allowed 
(Case  2). 


Year 

#  Open 
Water 
Days 

Mean  Annual 
Ice 

Thickness 

T(i)  Max 
°C 

T  Max 
°C 

H~p 

Diff 

cal/cm2/year 

Solar 

Net 

1 

33.7 

50.02 

-1.030 

-0.302 

2999.5 

260.1 

6745.8 

-8384.9 

1620.5 

2 

50  3 

45.19 

-0.329 

0.336 

2999.2 

315  7 

10131  0 

-11687.0 

1758  7 

3 

61.7 

39.95 

0.251 

0.862 

2998.9 

359.8 

12674.0 

- 14247.0 

1785.8 

4 

69.7 

35.94 

0.711 

1.287 

2998.7 

394.1 

14648.0 

-16254.0 

1786  7 

5 

77.3 

33.64 

1.027 

1  584 

2998.5 

420.0 

16062.0 

-18403.0 

1077.4 

6 

68.7 

37.30 

0.641 

1.221 

2998.4 

439.3 

14398.0 

-17006.0 

829  4 

7 

86.3 

33.71 

1.122 

1.676 

2998.4 

452.2 

16590.0 

-19625.0 

416.2 

8 

57.3 

43.91 

0.005 

0.646 

2998.3 

460.0 

11671.0 

-  15016.0 

114.1 

9 

74.7 

35.24 

0.866 

1.438 

2998.3 

463.7 

15382.0 

-17640.0 

1204.2 

10 

71  0 

36.12 

0.760 

1.331 

2998.3 

464.3 

14937.0 

-17488.0 

911.7 

11 

86.3 

32.50 

1.215 

1.768 

2998.3 

462.4 

17051.0 

-20072.0 

439  4 

12 

660 

39.22 

0.482 

1.078 

2998.3 

458.8 

13739.0 

- 16780.0 

416.1 

13 

86  7 

32  90 

1.151 

1.707 

2998.3 

453.7 

16769.0 

-19634.0 

587.9 

14 

67.7 

38.24 

0.574 

1.162 

2998.4 

447.6 

14135.0 

-17103.0 

478.2 

15 

87  7 

32.53 

1  184 

1.739 

2998  4 

440.8 

16915.0 

- 19928  0 

426.1 

16 

67.7 

38.47 

0.574 

1.163 

2998.5 

433.5 

14137.0 

-17136.0 

433  5 

17 

88.3 

32.37 

1.216 

1.770 

2998.5 

425.9 

17048  0 

-20087  0 

385  9 

18 

67.7 

38.44 

0.574 

1.163 

2998  6 

418.2 

14137.0 

-17187  0 

367.3 

19 

890 

32  24 

1.248 

1.802 

2998  6 

410.5 

17188.0 

-20241.0 

356.6 

20 

67  7 

38.45 

0.574 

1.164 

2998  6 

402.9 

14137.0 

-17229.0 

309  4 

T(1)  Max— maximum  temperature  at  depth  level  1  (2.5  m) 

T  Max— maximum  temperature  in  the  water  column 

H^— total  heat  input  into  the  water  column  by  advection  from  the  deep  ocean 
Diff— total  heat  input  into  the  water  column  by  diffusion  across  the  bottom 
Solar— total  heat  input  into  the  water  column  by  solar  radiation 
FB— total  heat  transfer  across  the  ice/ocean  or  ice/water  interface 
Net— +  Diff  +  Solar  +  F0 

fluxes  at  the  interfaces.  The  melt  at  the  surface  penetrating  solar  radiation  (Case  1G),  the  mean  annual 

increases,  since  solar  radiation  no  longer  penetrates  the  ice  thickness  was  320.47  cm  and  the  minimum  ice 

ice;  thus,  more  heat  is  available  for  melting  at  the  sur-  thickness  was  269.22  cm.  Case  2D  had  thicker  ice,  since 

face.  The  melt  at  the  bottom  increases,  since  some  of  some  of  the  heat  from  solar  radiation  was  stored  in 

the  heat  stored  in  the  water  column  is  released  and  the  heat  reservoir,  and  then  was  released  when  the  sur- 

melts  ice.  On  top  of  this  is  added  the  melt  due  to  the  face  temperature  falls  below  freezing.  For  this  case, 

heat  stored  in  the  heat  reservoir.  This  may  be  as  much  the  maximum  amount  of  heat  stored  in  the  reservoir 

as  that  melted  by  the  interface  heat  fluxes  alone.  Thus  during  the  first  year  was  1740.3  cal/cm5.  This  heat 

the  melting  of  ice  greatly  accelerates  until  the  heat  in  delayed  the  drop  in  the  surface  temperature  to  freez- 

the  reservoir  is  exhausted.  This  can  be  observed  from  ing  for  only  3  days,  but  greatly  affected  the  melt  of 

the  values  listed  in  Table  10.  This  process  is  repeated  ice  at  the  surface  during  the  ice-melt  phase.  The  total 

in  other  simulations  whenever  the  ice  thickness  becomes  amount  of  ice  melted  at  the  surface  was  66. 14  cm.  The 

so  thin  that  no  more  heat  can  be  stored  in  the  heat  total  growth  of  ice  at  the  bottom  was  42.49  cm.  The 
reservoir.  net  ice  melt  for  the  year,  then,  was  23.64  cm.  For 

When  the  ice  is  thicker,  the  heat  reservoir  will  Case  IG,  the  total  melt  at  the  surface  was  90.83  cm. 

increase  until  the  surface  temperature  of  the  ice  falls  The  total  growth  of  ice  at  the  bottom  was  44.03  cm, 

below  freezing.  To  examine  this  case,  the  initial  ice  and  the  net  ice  melt  for  the  year  was  46.81  cm.  The 

thickness  was  increased  to  340  cm  and  the  previous  amount  of  ice  grown  at  the  bottom  was  not  affected 

cases  were  repeated.  For  the  case  with  penetrating  solar  much  by  the  brine  pocket  simulation  in  this  case,  since 

radiation  (Case  2D)  the  mean  annual  ice  thickness  for  no  solar  radiation  penetrated  through  the  ice  to  warm 

year  1  was  331.42  cm,  and  the  minimum  ice  thickness  the  water  column.  The  main  difference  between  the 

that  year  was  293.35  cm.  For  the  case  without  cases  was  due  to  the  melt  at  the  surface. 


Table  9.  Listing  of  bottom  melt,  surface  melt,  total  melt,  and  ice  thickness  for  year  1  of  Case  1 F  and  Case  2 
from  the  onset  of  ice  melt  until  open  water  occurs  in  Case  2. 


Melt  at  Bottom 

Melt  at  Surface 

Total  Melt 

H 

Step  # 

C1F  C2 

C1F 

C2 

C1F  C2 

ClF 

C2 

502 

0.0703  0.0696 

0.0000 

0.0000 

0.0703  0.0696 

94.720 

91  723 

512 

0.0633  0.0771 

0.5712 

0.4681 

0.6346  0.5452 

88  363 

86  330 

522 

0.0649  0.0823 

0.5692 

0.4695 

0  6341  0.5518 

82.016 

80  836 

532 

0.0678  0.0874 

0.5607 

0.4649 

0.6285  0.5523 

75  701 

75.310 

542 

0.0715  0.0927 

0.5463 

0.4547 

0.6178  0.5474 

69  471 

69  809 

552 

0.0759  0.0984 

0.5263 

0.4394 

0.6022  0.5378 

63.375 

64  384 

562 

0.0812  0.1044 

0.5011 

0.4191 

0.5823  0.5235 

57  458 

59.081 

572 

0.0876  0.1109 

0.4712 

0.3944 

0  5588  0.5053 

51  761 

53.942 

582 

0.0950  0.1181 

0.4367 

0.3656 

0.5318  0.4837 

46  319 

49  005 

592 

0.1042  0.1251 

0.4028 

0.3370 

0.5071  0.4621 

41  143 

44  289 

602 

0.1151  0.1158 

0.3664 

0.4207 

0.4816  0.5365 

36.210 

36  226 

612 

0.1285  0  1395 

0.3250 

0.3607 

0  4535  0.5002 

31  547 

25  913 

622 

0.1450  0  1783 

0.2777 

0.2941 

0  4228  0  4724 

27  178 

20  128 

632 

0.1574  0.2189 

0.2318 

0.2175 

0.3892  0.4364 

23  133 

15  599 

642 

0  1812  0  2828 

0.1712 

0.1144 

0.3524  0.3972 

19  440 

1 1  448 

652 

0.2113  0  3619 

0.1012 

0.0000 

03125  0.3619 

16  133 

7  694 

662 

0.2501  0.3449 

0.0198 

0.0000 

0.2699  0  3449 

13  240 

4  164 

672 

0.2396  0.3211 

0.0000 

0.0000 

0  2396  0.3211 

10  721 

0  0000 

TOTAL  MELT: 

83  999 

91  723 

Table  10.  Listing  of  ice  melt,  heat  stored  in  the  heat  reservoir,  oceanic  heat  flux,  and  ice  thickness  for  year 
1  of  Case  2.  The  values  shown  occur  after  the  heat  in  the  reservoir  has  reached  its  maximum  and  is  being 
released. 

Step  # 

Melt  at  Bottom 

Melt  at  Surface 

Solar 

Total  Melt 

Fb 

H 

594 

0  1267 

0  3312 

640  35 

04579 

-  1  4807E-4 

43  37 

595 

0  1274 

0  3282 

644  60 

0  4556 

-  1  4827E-4 

42  91 

596 

0  1282 

0  3252 

625  39 

0  7786 

-  1  4847E-4 

42  14 

597 

0  1284 

0  4441 

593  41 

1  0166 

-  1  4590E-4 

41  12 

598 

0  1241 

0  4396 

561  76 

1  0034 

-  1  3117E-4 

40  12 

599 

0  1189 

0  4351 

530  43 

0  9891 

-  1  1379E-4 

39  13 

600 

0  1164 

0  4304 

499  44 

0  9771 

-  1  0307E-4 

38  15 

601 

0  1155 

0  4256 

468  80 

09668 

-9  6371E-5 

37  18 

602 

0  1158 

0  4207 

438  51 

0  9573 

-  9  21 12E-5 

36  23 

603 

0  1167 

04158 

400  17 

1  0646 

-  8  9333E-5 

35  16 

604 

0  1 1 79 

0  4104 

362  13 

1  0566 

-  8  6677E-5 

34  10 

605 

0  1194 

0  4049 

324  39 

1  0486 

-8  4113E-5 

33  06 

606 

0  1213 

0  3992 

286  91 

1  0410 

-8  241  IE-5 

32  01 

607 

0  1236 

0  3933 

249  69 

1  0339 

-8  1229E-5 

30  98 

608 

0  1263 

0  3872 

212  72 

1  0270 

-  8  0422E-5 

29  95 

609 

0  1291 

0  3810 

175  99 

1  0203 

-  7  9840E-5 

28  93 

610 

0  1323 

0  3745 

139  50 

1  0136 

-7  9451E-5 

27  92 

611 

0  1358 

0  3677 

103  25 

1  0069 

-7  9176E-5 

26  91 

612 

0  1395 

0  3607 

67  23 

1  0004 

-  7  8975E-5 

25  91 

613 

0  1435 

0  3534 

31  46 

0  9937 

-  7  8842E-5 

24  92 

614 

0  1478 

0  3457 

00 

0  9305 

-  7  8781E-5 

23  99 

615 

0  1532 

0  3379 

0  0 

0  4911 

-  8  1077E-5 

23  50 

616 

0  1574 

03320 

00 

04895 

-  8  5072E-5 

23  01 

617 

0  1612 

0  3260 

00 

0  4872 

-  8  7933E-5 

22  52 
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This  simulation  of  brine  pockets  is  rather  crude,  but 
it  shows  that  any  diversion  of  solar  heat  to  brine 
pockets  can  significantly  affect  the  growth  and  decay 
of  ice. 

7.  Case  3 — 

run  with  AIDJEX  data 

a.  Setup  of  model 

In  this  last  section,  the  ice/ocean  model  is  run  with 
forcing  data  obtained  from  the  AIDJEX  experiment. 
The  results  of  this  run  will  be  discussed  and  compared 
to  selected  fields  from  the  AIDJEX  data.  The  ice/ocean 
model  requires  the  following  as  input:  incoming  solar 
radiation,  incoming  longwave  radiation,  latent  heat 
flux,  sensible  heat  flux,  albedo,  winds,  snowfall  rates, 
and  the  heat  flux  from  the  deep  ocean.  The  model  also 
needs  the  initial  values  for  the  ice  thickness,  snow 
thickness,  and  the  temperature  and  salinity  profiles  in 
the  ocean.  Most  of  the  meterological  forcing  data  was 
obtained  from  two  AIDJEX  reports.  The  first  is  Report 
on  the  AIDJEX  Meterological  Experiment  by  Leavitt 
et  al.  (1978).  From  this  report,  the  latent  and  sensible 
heat  fluxes,  the  average  air  temperature,  and  the  wind 
speed  and  direction  were  obtained.  The  second  report 
is  Radiation  Program  during  AIDJEX:  A  Data  Report 
by  Hautzke  and  Hornof  (1978).  This  report  gives  the 
daily  total  of  the  incoming  solar  radiation,  and  the 
daily  average  tor  the  albedo.  The  incoming  longwave 
radiation  is  not  given  directly  in  the  AIDJEX  data,  but 
can  be  computed  from  the  average  air  temperature  with 
the  method  used  in  Parkinson  and  Washington  (1979). 
The  initial  ice  thickness  and  the  initial  temperature  and 
salinity  profiles  can  be  obtained  from  AIDJEX 
technical  reports.  Ice  thickness  measurements  were 
made  at  the  beginning  of  the  experiment.  These  values 
ranged  from  250  cm  to  470  cm  at  the  lour  camps:  the 
average  thickness  was  340  cm  This  value  is  used  to 
initialize  the  model.  The  snowfall  rates,  initial  snow 
thickness,  and  the  deep  ocean  heat  flux  are  not  given 
in  any  ol  the  AIDJEX  reports.  Therefore,  the  clima¬ 
tology  values  used  in  the  previous  model  runs  are  used 
m  the  runs  in  this  section. 

I  he  forcing  fields  from  the  AIDJEX  reports  are 
plotted  m  figures  35  through  41.  The  data  runs  from 
May  I.  1975.  to  April  29.  1976.  All  fluxes  except  for 
the  albedo  and  solar  radiation  are  given  at  6-hour  inter¬ 
vals.  Thus,  the  time  step  used  for  the  model  is  6  hours. 
The  solar  radiation  in  the  report  is  a  daily  total,  so 
the  question  arises  as  to  how  to  partition  the  solar 
radiation  data  into  6-hour  intervals  like  the  other  data 
Two  methods  of  partitioning  the  solar  radiation  will 
he  examined  The  first  method  is  to  assume  a  diurnal 
>.vdc  tor  the  solar  radiation.  There  is  no  readily 
available  data  on  the  diurnal  cycle  in  the  Beaufort  5>ca, 


but  one  can  be  deduced  from  Figure  10  in  the  report 
by  Pautzke  and  Hornof  (Fig.  42).  The  second  method 
is  to  evenly  distribute  the  daily  total  solar  radiation 
between  the  four  periods  within  a  day,  i.e.,  to  assume 
that  the  solar  flux  is  constant  throughout  the  day. 

The  heat  flux  from  the  deep  ocean  Hdeep  will  ini¬ 
tially  be  set  to  3.25  kcal/cm2/year;  however,  other 
values  will  be  tried.  A  listing  of  other  parameters  for 
this  simulation  is  shown  in  Table  11. 

b.  Diurnal  eye?? 

Most  of  the  input  data  is  given  at  0,  6,  12,  and  18 
hours  GMT.  Local  time  is  about  10  hours  earlier,  thus 
the  first  period  in  the  model  is  from  2  p.m.  to  8  p.m., 
the  second  period  is  from  8  p.m.  to  2  a.m.,  the  third 
period  is  from  2  a.m.  to  8  a.m.,  and  the  fourth  period 
is  from  8  a.m.  to  2  p.m.  The  first  and  fourth  periods 
are  daytime  hours,  and  the  second  and  third  periods 
are  nighttime  hours.  It  is  further  assumed  that  the  max¬ 
imum  solar  radiation  influx  occurs  during  the  fourth 
period,  i.e.,  the  8  a.m.  to  2  p.m.  interval.  Next,  the 
year  is  divided  into  four  groups.  The  first  group  is  sum¬ 
mer,  i.e,  the  months  of  May,  June,  and  July.  During 
these  months  some  solar  radiation  exists  at  each  period. 
In  particular,  30%  of  the  daily  total  of  solar  radia¬ 
tion  is  assumed  to  arrive  during  the  first  period,  10% 
during  each  of  the  next  two  periods  (the  nighttime 
hours),  and  50%  during  the  last  period.  The  second 


Table  1 1 .  Listing  of  parameters  and  assigned  values 
used  in  Case  3  tor  the  VML  model.  AIDJEX  data  is 
used  to  force  the  model. 


Variable  Mixed  Layer 

H 

340  cm 

30  cm 

F*> 

-  l  0306E-4  i  heal  change  (cal/cnv/sec) 

-  2  0°C 

T, 

Varies  with  salinity 

0 

l  3545E-12  cal/cm'1/  K‘/sec 

It 

064 

a  ^ 

0  07 

q. 

72  cal/cm' 

q. 

72  cal/cm  ‘ 

f 

i  0 

MID 

Variable 

SFF 

3  341E  S  cm/sec 

Time  Step 

6  hours/30  minutes 

Forcing 

Heal  Fluxes 

AIDJEX  Data 

Winds 

AIDJEX  Data 

T4S  Protiles 

AIDJEX  Data 
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group  is  the  fall  and  spring  months  of  August, 
September,  March,  and  April.  During  these  months 
no  solar  radiation  arrives  during  the  nighttime  hours. 
The  total  solar  radiation  for  the  day  is  divided  between 
the  daytime  hours:  30%  for  one,  70%  for  the  other. 
The  third  group  is  the  months  of  October  and 
February.  For  these  months,  solar  radiation  exists  dur¬ 
ing  only  one  of  the  daytime  periods.  The  rest  of  the 
time  the  solar  radiation  is  zero.  The  last  group  is  the 
winter  months  of  November,  December,  and  January, 
where  the  solar  radiation  is  zero  throughout  the  day. 
This  partitioning  of  the  solar  radiation  is  shown 
graphically  in  Figure  43. 

The  time  series  for  the  solar  radiation  field,  which 
results  when  the  above  partitioning  is  applied  to  the 
solar  radiation  data  in  Figure  35  is  shown  in  Figure  44. 
The  resulting  solar  radiation  values  are  quite  variable, 
with  most  of  the  energy  being  concentrated  during  the 
daytime  hours  and  little  or  none  existing  during  the 
nighttime  hours. 

c.  No  diurnal  cycle 

The  second  method  of  partitioning  the  solar  radia¬ 
tion  is  to  assume  that  the  solar  flux  during  the  day  is 
constant.  The  solar  flux  is  then  determined  by  dividing 
the  daily  total  solar  radiation  by  the  seconds  in  a  day. 
The  resulting  time  series  for  this  partitioning  of  the 
solar  radiation  is  shown  in  Figure  45.  The  extremes 
in  the  values  are  much  less  in  this  case,  and  the  energy 
is  more  evenly  distributed  between  time  steps  than  in 
the  diurnal  case. 

d.  Diurnal  cycle  versus  nondiurnal  cycle 

To  examine  how  the  partitioning  of  the  solar  radia¬ 
tion  may  affect  the  model  results,  the  ice/ocean  model 
was  run  for  1  year  with  the  diurnal  cycle  and  for  1  year 
without  the  diurnal  cycle.  The  model  results  for  these 
cases  will  be  discussed  in  this  section. 

Table  12  is  a  listing  of  some  results  from  the  model 
simulations.  Tne  case  with  the  diurnal  cycle  has 
somewhat  thinner  ice  than  the  case  without.  This  dif¬ 
ference  in  the  ice  thickness  is  due  to  increased  melting 
at  the  surface.  Growth  of  ice  at  the  bottom  is  somewhat 
larger  in  the  diurnal  case,  but  not  enough  to  offset  the 
ice  melt  at  the  surface.  The  net  result  is  increased 
melting  and  thinner  ice. 

The  surface  temperatures  throughout  the  year  for 
both  simulations  are  shown  in  Figures  46  and  47.  In 
these  figures,  it  is  obvious  that  the  temperature 
extremes  are  greater  in  the  diurnal  case.  There  can  be 
as  much  as  30°C  difference  in  the  daytime  and  the 
nighttime  temperatures  in  this  case  (this  value  occurred 
in  April)  The  maximum  difference  in  the  nondiurnal 
case  was  about  8°C.  For  the  diurnal  case,  Figure  46, 
Ts  first  reaches  zero  degrees  on  May  5,  and  fre¬ 
quently  gets  this  warm  throughout  May  and  the  rest 
of  the  summer  (in  the  model  if  7\  is  above  the 


Table  12.  Listing  of  results  from  a  1-year  simulation 
of  the  VML  model  forced  with  data  from  AIDJEX. 


No  Diurnal  Cycle 

Diurnal  Cycle 

Mean  annua! 
ice  thickness 
(cm) 

325  99 

310.01 

Mean  annual 
oceanic  heat  flux 
(cat/cm2/sec) 

- 1 .1 146E-4 

-1  1065E-4 

Deep  ocean  heat 
input  (cal/cm2/year) 

3291.90 

3291  90 

Diffusive  heat  input 
(cal/cm2/year) 

264  1  2 

264  12 

Solar  heat  input 
(cal/cm2/year) 

27  05 

40  78 

Total  oceanic  heat 
flux  (cal/cm2/year) 

3514  90 

3489  40 

Net  heating 
(cal/cm2/year) 

68  13 

107  37 

Number  of  snow- 
free  days 

63  25 

76  00 

Net  change  in  ice 
thickness  for  the 

year  (cm) 

2  32 

-  17  07 

Ice  melt  at  surface 
for  the  year  (cm) 

16  77 

42  33 

Ice  growth  at  bottom 
for  the  year  (cm) 

1909 

25  25 

freezing  point,  it  is  set  to  the  freezing  temperature,  and 
the  additional  heat  is  applied  to  ice  melt).  However, 
the  surface  temperature  for  the  nondiurnal  case,  Fig¬ 
ure  47,  does  not  reach  0°  until  the  end  of  May.  Thus, 
surface  melting  begins  much  earlier  in  the  diurnal 
simulation  than  in  the  nondiurnal  simulation.  In  fact, 
the  snow  cover  is  totally  melted  by  step  189  (June  16) 
in  the  diurnal  case,  but  not  until  step  240  (June  29) 
in  the  nondiurnal  case.  Thus,  ice  melt  at  the  surface 
begins  about  13  days  earlier  in  the  diurnal  case.  From 
the  figures,  it  is  also  seen  that  during  June  and  July 
the  surface  temperature  for  the  nondiurnal  case 
remains  at  0°  most  of  the  time,  but  that  the  tempera¬ 
tures  for  the  diurnal  case  decrease  at  night.  For  the 
nondiurnal  case  the  melting  at  the  surface  will  be  fairly 
constant,  while  for  the  diurnal  case,  melting  will  be 
strong  during  the  day  and  weak  at  night.  In  fact,  on 
some  nights  ice  forms  on  the  bottom  in  the  diurnal  case. 

Ice  melt  at  the  surface  of  the  ice,  and  ice  melt  and 
growth  at  the  bottom  for  the  two  simulations  are  shown 
in  Figures  48  through  51.  Ice  melt  at  the  surface  in  both 
cases  lasts  fro.n  June  through  August.  The  magnitude 
of  ice  melt  at  the  surface  is  generally  largest  during 
the  month  of  July.  The  maximum  melt  per  time  period 
in  the  diurnal  case  is  just  over  1  cm,  while  for  the  non 
diurnal  case  it  is  less  than  0.5  cm.  The  figures  of  growth 
or  melt  of  ice  at  the  bottom  (Figs.  50  and  51)  show 
the  expected  trend  of  melt  during  the  summer  and 


growth  during  the  winter  in  both  cases.  For  the  diur¬ 
nal  case,  however,  melt  may  occur  during  the  daytime 
hours,  while  growth  occurs  during  the  nighttime  hours. 
During  the  spring  and  fall  the  melt  during  the  day  can 
almost  be  offset  by  growth  at  night.  But,  as  mentioned 
earlier,  the  surface  melt  in  the  diurnal  case  is  sufficient 
to  result  in  thinner  ice  in  the  diurnal  case  than  in  the 
nondiurnal  case. 

Another  factor  that  affects  ice  melt  at  the  surface 
is  the  brine  pockets  simulation  discussed  in  the  previous 
section.  The  storage  of  heat  in  the  heat  reservoir  occurs 
during  the  summer  when  the  ice  is  snow-free.  In  the 
diurnal  case,  Ts  generally  falls  below  freezing  at  night 
during  the  summer.  Thus,  no  heat  is  diverted  into  the 
heat  reservoir  during  this  time,  instead  most  if  not  all 
of  the  heat  stored  in  the  reservoir  during  the  day  will 
be  released  at  night.  The  total  amount  of  heat  stored 
in  the  reservoir  during  a  day’s  time,  then,  is  less  in  the 
diurnal  case  than  in  the  nondiurnal  case.  The  maximum 
amount  of  heat  stored  in  the  reservoir  at  any  one  time 
is  about  30  cal/cm:  in  the  diurnal  case.  While  for  the 
nondiurnal  case,  the  heat  in  the  reservoir  steadily  in¬ 
creases  to  a  maximum  of  about  500  cal/cm2,  and 
then  decreases  to  zero  when  the  surface  temperature 
stays  below  freezing  (Figs.  52  and  53).  Thus,  the  simu¬ 
lation  of  brine  pockets  has  less  effect  on  the  ice  thick¬ 
ness  when  a  diurnal  cycle  is  assumed  for  the  solar  flux 
than  when  the  solar  flux  is  constant  throughout  the 
day.  More  information  on  brine  pockets  is  needed  to 
determine  the  importance  of  this  effect  and  how  best 
to  simulate  it. 

e.  Model  results  versus  AIDJEX  data 

The  results  of  a  1  year  simulation  of  the  ice/ocean 
model  forced  with  AIDJEX  forcing  data,  as  described 
in  Section  "'a.  is  compared  with  data  obtained  during 
AIDJEX  in  this  section.  The  diurnal  cycle  for  the  solar 
radiation  is  chosen  for  this  run  (Case  3),  and  the  ini¬ 
tial  temperature  and  salinity  profiles  are  from  May  10, 
19"'5,  Station  Blue  Fox  (Figs.  54  and  55). 

Comparison  will  be  made  on  how  well  the  model 
simulates  the  mixed-layer  characteristics  rather  than 
on  the  growth  and  decay  of  ice.  Considerable  data  exist 
as  to  the  temperature  and  salinity  structure  in  the  water 
column  below  the  ice  during  AIDJEX,  but  very  little 
on  the  ice  thickness,  and  ice  growth  and  decay. 

The  mixed-layer  temperature  computed  by  the  model 
for  I -year  is  shown  in  Figure  56.  For  comparison,  the 
mixed-layer  temperature  at  Station  Blue  Fox  and  Sta¬ 
tion  Snowbird  are  shown  in  Figures  57  and  58.  From 
these  figures,  one  sees  that  the  ice/ocean  model 
simulates  the  general  trends  of  temperature  increase 
in  the  summer  and  decrease  in  the  fall  and  winter  rather 
well  The  magnitude  of  the  mixed-layer  temperature 
is  also  reasonable.  The  extremes  in  the  temperature, 
however,  are  not  well  simulated.  The  data  indicates 
that  the  mixed-layer  temperature  is  warmer  in  the 


summer  than  the  model  indicates  and  is  cooler  in  the 
winter.  Similar  plots  for  the  mixed  layer  salinity  are 
shown  in  Figures  59  through  61.  The  model  does  a 
somewhat  better  job  simulating  the  salinity,  except  that 
the  salinity  is  a  little  too  fresh  in  the  winter  and  early 
spring.  Finally,  the  mixed-layer  depth  determined  from 
the  salinity  (as  described  earlier  in  this  paper)  is  plot¬ 
ted  in  Figures  62  through  64.  The  model  again  follows 
the  general  trends  quite  well.  It  causes  shallowing  of 
the  mixed-layer  in  the  summer  and  deepening  in  the 
winter.  The  model,  however,  shallows  too  much  in  the 
summer  and  underestimated  the  mixed-layer  deepen¬ 
ing  in  the  winter.  This  discrepancy  may  be  due  to 
several  factors.  The  turbulence  parameters  prescribed 
in  the  model  may  need  to  be  tuned  to  cause  the  mixed 
layer  to  deepen  more.  The  shallow  mixed  layer  may 
also  be  due  to  uncertainties  in  some  of  the  forcing 
values  used  in  the  model,  for  example,  the  snowfall 
rate  and  the  heat  flux  from  the  deep  ocean.  The  dif¬ 
ferences  may  also  be  due  to  advection,  which  is  not 
accounted  for  in  the  model.  Some  of  these  possibilities 
will  be  examined  later  in  this  report. 

A  comparison  of  the  temperature  and  salinity  pro¬ 
files  for  the  model  results.  Station  Blue  Fox,  and  Sta¬ 
tion  Snowbird  is  shown  in  Figures  65  through  76.  The 
solid  line  in  the  plots  is  the  model-determined  profile, 
and  the  dashed  line  is  from  AIDJEX  data.  Again  the 
general  trends  and  magnitudes  compare  well,  except 
that  the  mixed  layer  is  too  shallow  and  the  profiles  from 
the  data  show  more  variability.  In  the  summer,  the 
upper  ocean  (30  m  to  50  m)  is  warmer  in  the  tempera¬ 
ture  profiles  from  the  data  than  in  the  model  profiles. 
In  the  winter,  the  temperature  of  the  upper  ocean  is 
somewhat  cooler  in  the  data  profiles  than  in  the  model 
profiles.  The  salinity  profiles  compare  better,  but  again 
the  mixed  layer  is  too  shallow,  and  the  salinity  tends 
to  be  a  little  fresher  in  the  model  profiles  than  in  the 
data  profiles,  particularly  when  compared  to  Station 
Blue  Fox.  Maximums  and  minimums  for  the  mixed- 
layer  temperature,  salinity,  and  mixed-layer  depth  for 
the  four  AIDJEX  camps  and  for  several  model  runs 
are  shown  in  Table  13. 

In  the  depth  range  of  50  m  to  100  m,  the  temperature 
profiles  from  the  model  are  in  rough  agreement  with 
data,  but  the  model  does  not  indicate  the  temperature 
maximum  around  60  m  depth,  which  is  quite  apparent 
at  both  Blue  Fox  and  Snowbird  in  the  spring  and  sum¬ 
mer.  The  model  also  does  not  simulate  well  the 
variability  in  the  temperature  in  this  depth  range. 

The  temperature  maximum  that  occurs  around  60  m 
depth  in  the  temperature  profiles  from  data  is  due  to 
Pacific  water,  which  enters  the  Arctic  through  the 
Bering  Strait  and  eventually  flows  into  the  Beaufort 
Sea.  This  Pacific  water  shows  up  in  the  profiles  from 
50  m  to  130  m  depth.  The  temperature  and  salinity  pro¬ 
files  from  AIDJEX  data  are  also  perturbed  by  eddies. 
Several  eddies  were  observed  during  the  AIDJEX 
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Table  13.  Listing  of  maximums  and  minimums  in  temperature,  salinity,  and 
mixed-layer  depth  for  the  four  camps  of  AIDJEX  and  four  VML  model  runs. 


Mixed-Layer  Characteristics 

Bluetox  Snowbird  Big  Bear  Caribou 


Max 

Min 

Max 

Min 

Max 

Min 

Max 

Min 

Temperature 

-1.20 

-  1.73 

-1.44 

-  1.78 

-1.46 

-  1.71 

-  1.40 

-  1  76 

Salinity 

31.32 

29.77 

30.98 

29.59 

31.17 

29  76 

30  48 

29  79 

Depth 

65.0 

10.0 

60.0 

10.0 

60.0 

15.0 

60.0 

15.0 

Model— 

-Case  3 

Model— 

Case  3J 

Model- 

Case  30 

Model— 

Case  30 

Max 

Min 

Max 

Min 

Max 

Min 

Max 

Min 

Temperature 

- 1 .616 

-1.660 

-1  614 

-1  665 

-1.621 

-1  668 

-1.013 

-1  689 

Salinity 

30.598 

29.794 

30.674 

29  761 

30.741 

29  883 

31.114 

29.460 

Depth 

50.0 

2.5 

50  0 

2  5 

50.0 

25 

500 

2.5 

experiment.  These  eddies  were  located  from  50  m  to 
300  m  depth.  Thus,  advection  is  important  in  simu¬ 
lating  the  temperature  and  salinity  structure  in  the 
Beaufort  Sea,  and  probably  accounts  for  some  of  the 
differences  between  model  results  and  data. 

As  mentioned  earlier,  uncertainties  in  some  of  the 
forcing  fluxes  may  also  account  for  some  of  the 
discrepancies  in  the  temperature  and  salinity  fields. 
Two  easily  examined  fluxes  are  the  snowfall  rate  and 
the  heat  flux  from  the  deep  ocean.  To  test  how  the 
snow  cover  might  affect  the  ocean  under  the  ice,  the 
snowfall  rate  was  divided  by  two,  and  the  initial  snow 
thickness  value  was  decreased  to  20  cm.  Snow  insulates 
the  ice  and  the  ocean  from  the  atmosphere;  thus,  less 
snow  means  less  of  an  insulating  effect.  In  this  test  case 
(Case  3J),  the  result  was  an  increase  in  the  growth  of 
ice  in  the  winter  and  less  melting  in  the  summer.  About 
10  cm  more  ice  grew  in  the  winter  in  this  case  than  in 
the  previous  case.  This  had  a  small,  but  noticeable 
effect  on  the  mixed-layer  depth.  The  mixed  layer 
deepened  to  32  m  instead  of  25  m  (Fig.  77).  Similar 
results  occurred  when  the  heat  flux  from  the  deep  ocean 
is  reduced  to  1.5  kcal/cm;/year  (Case  3Q;  Fig.  78). 
Both  changes  increased  the  growth  of  ice  in  the  winter, 
which  results  in  increased  salinity  flux  at  the  surface 
and  increased  deepening  of  the  mixed  layer. 

f.  Open  water  and  AIDJEX 

The  AIDJEX  forcing  and  initial  conditions  do  not 
yield  open  water  during  the  simulation.  This  is  not  sur¬ 
prising,  since  the  measurements  were  taken  on  ice 
flows.  No  direct  information  about  open  water  dur¬ 
ing  AIDJEX  is  available,  except  that  camp  Big  Bear 
broke  up  and  had  to  be  abandoned.  It  would  be  inter¬ 
esting,  however,  to  see  how  an  open-water  period  might 
affect  the  temperature  and  salinity  structure  in  the 
water  column  using  the  same  forcing  as  in  the  previous 
AIDJEX  cases. 


To  force  open  water,  the  ice  thickness  was  reduced 
to  65  cm,  and  the  heat  flux  from  the  deep  ocean  Hdefp 
was  increased  to  6.0  kcal/cm-Vyear  (Case  3N).  The 
result  was  an  open-water  period  of  4.25  days.  Another 
test  was  run  with  an  initial  ice  thickness  of  63  cm 
(Case  30).  The  result  was  an  open-water  period  of 
43  days.  The  mixed-layer  temperature,  mixed-layer 
salinity,  and  the  ice  thickness  for  both  cases  are  shown 
in  Figures  79  through  81.  These  cases  yield  similar 
results  up  to  the  open-water  period.  During  the  open- 
water  period,  the  temperature  for  Case  30  (the  solid 
line  in  the  plots)  warms  to  -1.013  C,  while  the 
temperature  for  Case  3N  (the  dashed  line  in  the  plots) 
warms  to  -  1 ,598°C.  The  salinity  for  Case  30  is  fresher 
than  for  Case  3N,  during  and  following  the  open-water 
period.  After  the  open-water  period,  the  temperature 
and  salinity  values  for  both  cases  gradually  approach 
the  same  value.  The  ice  thickness  for  Case  3N  increases 
rapidly  after  the  open-water  period,  and  by  the  end 
of  the  simulation  is  greater  than  the  ice  thickness  in 
Case  30.  This  rapid  increase  in  ice  thickness  is  due  to 
thinner  snow  cover.  This  effect  of  snow  cover  was 
discussed  in  Section  3b.  The  salinity  mixed-layer  depth 
for  Case  30  is  shown  in  Figure  82.  The  mixed-layer 
depth  plot  for  Case  3N  is  similar,  except  that  the  mixed 
layer  starts  to  deepen  sooner,  around  May  27.  The 
mixed-layer  depth  in  both  cases  reaches  50  m  in  the 
winter.  The  rapid  growth  of  ice  following  the  open- 
water  period  dumps  large  amounts  of  salt  and  causes 
the  mixed  layer  to  deepen. 

The  temperature  and  salinity  profiles  for  Case  30 
compared  with  Blue  Fox  are  shown  in  Figures  83 
through  88.  In  this  case,  the  temperature  of  the  upper 
50  cm  is  closer  to  the  temperature  in  Blue  Fox  than 
was  true  in  Case  3,  at  least  for  most  of  the  year.  In 
late  August,  open  water  occurs,  and  by  September  1 
the  surface  waters  have  warmed  to  about  -  1 .0°C.  A 
warm  spot,  such  as  the  one  observed  in  the  climatology 
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case,  is  formed  in  this  case.  The  warm  spot  decreases 
as  the  mixed  layer  deepens,  and  by  late  December  the 
profiles  for  Case  30  and  Blue  Fox  look  very  much 
alike.  The  salinity  profiles  compare  well  also,  mainly 
after  the  open-water  period. 

It  is  interesting  that  the  open-water  period  creates 
a  temperature  maximum  of  about  the  same  magnitude 
and  at  about  the  same  depth  as  the  temperature  max¬ 
imum  from  the  Pacific  water  This  is  probably  a  coinci¬ 
dence,  but  only  more  research  on  the  relationship 
between  open  water  and  ice  growth  and  movement  will 
tell  whether  open  water  accounts  for  any  of  the  features 
observed  in  the  AIDJEX  profiles. 

8.  Summary  and  conclusions 

A  one-dimensional  sea-ice  and  ocean  mixed-layer 
model  has  been  developed  and  was  used  to  investigate 
the  interaction  between  ice  and  ocean.  This  model, 
forced  with  climatology  forcing  fields,  was  compared 
to  Semtner’s  O-layer  ice  model.  With  the  same  forcing 
and  initial  conditions,  Semtner’s  model  yielded  a  6-year 
cycle  of  open  wate  .  Our  model  yielded  open  water 
nearly  every  year.  The  primary  difference  between  the 
models  is  the  treatment  of  the  mixed  layer  in  the  ocean. 
In  our  ice/ocean  model,  the  depth  of  the  mixed  layer 
varies  in  response  to  forcing.  The  exchange  of  heat 
between  the  ocean  and  the  ice  varies  with  the  mixed- 
layer  depth  and  heat  input  from  below.  The  heat 
exchange  was  the  greatest  when  the  mixed  layer  deep¬ 
ened  and  entrained  warmer  water  from  below  the 
mixed  layer.  Warmer  water  existed  below  the  mixed 
layer  because  the  temperature  tended  to  increase  with 
depth  in  the  simulation.  Large  temperature  differences 
in  the  water  column  developed  when  several  days  of 
open  water  occurred.  Solar  radiation  penetrating  into 
the  ocean  would  warm  the  subsurface  layers,  forming 
a  warm  spot  in  the  water  column.  Some  heat  remained 
in  this  warm  spot  even  after  ice  was  re-formed.  This 
occurred,  because  the  water  column  was  highly  strati¬ 
fied  during  this  time,  and  the  exchange  of  heat  between 
layers  was  slow .  The  heat  in  the  warm  spot  dissipated 
as  the  ice  grew,  dumping  salt  and  deepening  the  mixed 
layer.  This  heat  was  often  released  in  large  bursts, 
which  cither  slowed  the  freezing  rate  or  melted  a  few 
centimeters  of  ice.  The  rate  at  which  the  heat  was 
released  trom  the  mixed  layer  was  tied  to  the  rate  of 
mechanical  stirring  generated  by  shear  stresses.  This 
simulation  of  the  ice/ocean  model  showed  that  a 
variable  depth  mixed  layer  can  significantly  alter  the 
results  of  a  sea-ice  model.  This  ice/ocean  model  should 
improve  sea-ice  forecasting  by  better  simulating  the 
mixed  layer.  The  biggest  improvement  would  be  in  the 
marginal  ice  zones  where  open  water  frequently  occurs. 

The  sensitivity  of  the  ice/ocean  model  to  forcing  and 
various  parameters  was  also  investigated.  Changes  in 
the  snow  cover,  assumed  flux  of  heat  from  the  deep 


ocean,  and  incoming  longwave  radiation  were  shown 
to  alter  the  growth  of  ice,  which  in  turn  increased  or 
decreased  the  deepening  of  the  mixed  layer  and  the 
exchange  of  heat  between  the  ice  and  ocean.  Large 
effects  occurred,  depending  on  the  treatment  of  the 
solar  flux.  The  brine  pockets  simulation  showed  that 
diversion  of  solar  radiation  into  the  brine  pockets  could 
significantly  decrease  the  melt  at  the  surface  of  the 
snow /ice.  This  simulation,  however,  resulted  in  rapid 
ice  melt  if  the  ice  was  thin,  and  the  maximum-allowed 
storage  of  heat  in  the  brine  pockets  was  reached.  Some 
solar  radiation  was  allowed  to  penetrate  the  ice  and 
warm  the  water  column  below  the  ice.  This  increased 
the  melt  of  ice  at  the  bottom,  but  the  surface  melt 
dominated.  The  brine  pockets  effect  was  large  in  the 
first  simulation,  since  the  solar  flux  was  constant 
throughout  the  day.  The  surface  temperature  did  not 
drop  below  freezing  during  most  of  the  summer.  Thus, 
the  heat  in  the  reservoir  increased  each  time  step  and 
diverted  heat  from  surface  melt.  If  a  diurnal  cycle  was 
assumed  for  the  incoming  solar  radiation,  the  brine 
pockets  simulation  had  little  effect  on  ice  melt.  The 
surface  temperature  frequently  dropped  below  freez¬ 
ing  at  night  during  the  summer.  Thus,  no  heat  was 
diverted  into  the  reservoir,  and  any  heat  in  the  reser¬ 
voir  was  released  at  night.  The  diurnal  cycle  also  con¬ 
centrated  more  heat  during  the  day,  resulting  in 
increased  surface  melt  and  frequent  ice  growth  at  night. 
Surface  melt  during  the  summer  was  sufficient  to  yield 
thinner  ice  than  when  no  diurnal  cycle  was  assumed. 

The  ice/ocean  model  was  forced  with  data  obtained 
during  the  AIDJEX  experiment.  The  results  of  this 
simulation  showed  that  the  ice/ocean  model  did  a  good 
job  of  simulating  the  general  trends  in  the  mixed-layer 
characteristics.  The  model,  however,  underpredicted 
the  deepening  of  the  mixed  laver  and  missed  the  high 
variability  observed  in  the  ocean  in  the  Beaufort  Sea. 
The  omission  of  advection  from  the  model  is  a  primary 
reason  for  many  of  the  discrepancies,  since  the 
temperature  and  salinity  fields  in  the  Beaufort  Sea  are 
perturbed  by  Pacific  water  and  transit  eddies.  Other 
possibilities  for  the  discrepancies  were  uncertainties  in 
some  of  the  forcing  fields,  such  as  the  snowfall  and 
the  heat  input  from  the  deep  ocean.  These  forcing  fields 
affect  the  growth  of  ice,  which  in  turn  affects  the 
mixed-layer  response. 

An  open-water  simulation  was  examined  for  the 
AIDJEX  case.  The  open-water  period  lasted  for 
43  days.  The  surface  temperature  of  the  water  column 
warmed  to  -  1.013  C  during  this  time,  and  a  warm 
spot  formed  around  10  m  depth.  This  warm  spot  cooled 
and  migrated  down  the  water  column  as  ice  re-formed 
and  grew.  By  the  end  of  the  year  a  warm  spot  similar 
to  that  attributed  to  the  Pacific  water  was  formed  in 
the  model.  This  result  was  interesting,  but  more 
research  is  required  before  we  can  put  much 
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significance  to  it,  since  we  have  no  data  about  open 
water  in  the  area  during  AIDJEX. 

This  ice/ocean  model  has  shown  that  a  variable 
depth  mixed-layer  ocean  model  can  simulate  the  mixed 
layer  in  the  Arctic,  but  the  model  should  be  improved 
by  including  advection.  In  the  future,  we  also  plan  to 
develop  a  three-dimensional  version  cf  this  model  and 
to  improve  the  simulation  of  the  ocean  currents  and 
shear  stresses. 
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Figure  1.  Vertical  grid  for  variable  mixed  layer  ice /ocean  model.  Ts  is  the  temperature  at  the  surface  of 
the  snow/ice  system.  TB  is  the  temperature  at  the  bottom  of  the  ice.  H|  and  Hs  are  the  ice  and  snow 
thicknesses,  respectively.  The  quantities  T,  S,  u,  and  v  for  the  ocean  are  defined  at  the  depths  indicated 
in  the  figure.  All  turbulence  quantities  are  defined  midway  between  these  depths. 
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Figure  2.  Temperature  profile  from  Station  Blue  Fox,  January  1,  1976.  This  profile  is  used  to  initialize 
the  ice/ocean  model. 
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Figure  4.  Solar  radiation  fluxes  used  to  force  the  climatology  test  cases.  Fluxes  are  determined  from  Fletcher’s 
(1965)  climatology  for  the  Arctic. 
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Figure  7.  Sensible  heat  fluxes  used  to  force  the  climatology  test  cases.  Fluxes  are  determined  from  Fletcher’s 
(1965)  climatology  for  the  Arctic. 
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Figure  8.  V -component  of  the  wind  velocity  determined  from  monthly  averaged  winds  from  FNOC’s  analysis. 
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Figure  9.  V -component  of  the  wind  velocity  determined  from  monthly  averaged  winds  from  FNOC’s  analysis. 
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Figure  12.  The  yearly  variability  of  the  ice  and  ice  plus  snow  thickness  for  a  20-year  simulation  of  the  CML 
model  and  the  VML  model.  The  upper  plot  is  the  results  from  the  CML  model.  The  lower  plot  is  the  results 
from  the  VML  model. 
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Figure  13.  The  mixed-layer  temperature  for  year  2  from  the  CML  simulation. 
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Figure  14.  The  mixed-layer  temperature  for  year  2  from  the  VML  simulation. 
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f  igure  15.  The  mixed-layer  depth  for  year  2  from  the  VML  simulation.  The  mixed-layer  depth  is  the 
imum  depth  at  which  the  salinity  has  changed  by  no  more  than  0.1  ppt.  from  its  surface  value. 
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Figure  16.  Some  temperature  prof ile  plots  from  the  VML  simulation.  These  profiles  were  taken  just  before , 
during,  and  just  after  the  open-water  period  in  year  2  of  the  simulation. 
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Figure  19.  Some  temperature  profile  plots  from  the  VML  simulation.  These  profiles  were  taken  just  before, 
during,  and  just  after  the  open-water  period  in  year  10  of  the  simulation.  Open  water  exists  from  step  659 
to  step  780. 
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Figure  20.  Some  temperature  profile  plots  from  year  10  of  the  VML  simulation.  These  prof iles  were  taken 
following  the  open  water  period,  during  the  time  ice  is  growing. 
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Figure  21.  Some  temperature  profile  plots  from  year  1 1  of  the  VML  simulation.  These  profiles  were  taken 
from  the  first  half  of  the  year  from  January  25  to  July  20. 
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Figure  22.  The  oceanic  heat  flux  from  the  VML  simulation  for  year  10. 
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Figure  23.  The  mixed-layer  depth  determined  from  the  salinity  for  year  10  of  the  VML  simulation. 
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Figure  24.  Plots  for  the  oceanic  heat  flux,  the  salinity  mixed-layer  depth,  and  the  ice  thickness  for  the  last 
5  years  of  the  VML  simulation. 
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flux  for  year  10  of  the  VML  simulation  for  the  case  when  the  winds  were  set 
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Figure  27.  Temperature  profiles  of  the  values  used  in  the  calculation  of  sound  speed . 
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Figure  2%.  Salinity  profiles  of  the  values  used  in  the  calculation  of  sound  speed. 
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/  inure  29  Ice  thickness  and  oceanic  heat  flux  values  for  a  20-year  simulation  of  the  VML.  model  when 
H  /  X  kcai  cm-  year. 
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Figure  31.  Ice  thickness  and  oceanic  heat  flux  values  for  a  20-year  simulation  of  the  VML  model  when 
Hdcep  =  4.0  kcal/cm2 /year. 
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Figure  32.  Ice  thickness  and  oceanic  heat  flux  values  for  a  20-year  simulation  of  the  VML  model  when 
H<jeep  -6-0  kcal/cm2/year. 
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Figure  33.  ke  thickness  and  oceanic  heat  flux  values  for  a  20-year  simulation  of  the  VML  model  with  brine 
pockets  included  (Case  2). 
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Figure  34.  Comparison  of  temperature  prof iles  for  Case  IF  and  Case  2.  Case  IF  is  without  the  brine  pockets 
simulation,  and  Case  2  is  with  the  brine  pockets. 


TEMPERATURE  PROFILE  PLOTS  TOR  MODEL  RUN  CASES  IP  AND  2 
TERR  NO.  1 


CRl/CM*  *2/DRY 


800.00 


SOLAR  ROD  I  AT  I  ON 
AIDJEX 


CPL/CM*  *2/SEC 


LONGWAVE  RADIATION 
AIDJEX 


I.C  146.9  292.8  433.7  584.6  "30.5  076.4  1022.3  1168.2  1314.1  1460.0 

5  01 ' 75  9/01/75  1/01/76  5/0!  -6 

Figure  37.  Longwave  radiation  fluxes  from  May  1,  1975,  to  April  29,  1976.  These  values  were  computed 
from  the  average  air  temperature  obtained  during  AIDJEX,  using  the  method  shown  in  Parkinson  and 
Washington  (1979). 


58 


62 


Figure  42.  Three-dimensional  perspective  of  the  flux  of  short-wave  radiation  over  one  year.  < From  Pautzke 
and  Hornof,  1978 ) 
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Figure  44.  The  solar  radiation  field  which  results  when  the  diurnal  cycle  shown  in  Figure  43  is  applied  to 
the  daily  total  solar  radiation  shown  in  Figure  35. 
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Figure  45.  The  solar  radiation  field  that  results  when  the  solar  flux  is  assumed  constant  throughout  the 
day  and  applied  to  the  daily  total  solar  radiation  shown  in  Figure  35. 


66 


TEMPERATURE  fC) 


diurnal  cycle 


SURERCE  TEMPERATURE 


SURFACE  TEMPERATURE 


SURFACE  TEMPERATURE 


215  257  269  261  293  306  317  329  331  353  365 

JHN  TCB  MRR  RPR 


Figure  46.  A  time  series  of  the  surface  temperature  computed  by  the  VML  model,  when  the  diurnal  cycle 
is  assumed  for  the  solar  flux. 
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f  igure  47.  A  time  series  of  the  surface  temperature  computed  by  the  VML  model  when  the  constant  daily 
flux  is  assumed  for  the  solar  flux. 


68 


01URNRL  CYCLE 


ICE  MELT  AT  SURFACE 


:rc  OCT  NOV  DCC 


ICC  MELT  AT  SURrACC 


Figure  48.  Ice  melt  at  the  surface  of  the  ice  when  the  diurnal  cycle  is  assumed. 
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Figure  49.  Ice  melt  at  the  surface  of  the  ice  when  no  diurnal  cycle  is  assumed. 
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Figure  51.  Ice  growth  or  melt  at  the  bottom  of  the  ice  when  no  diurnal  cycle  is  assumed. 
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/  inure  53.  Heat  stored  in  the  brine  pocket  heat  reservoir  when  no  diurnal  cycle  is  assumed. 
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Figure  54.  Temperature  profile  from  May  10,  1975,  camp  Blue  Fox.  This  profile  is  used  to  initialize  the 
VML  model. 


JErTH 


MODEL  PROFILES 
SALINITY  STEP  1 


09.5  30.0  30.5  31.  C  31.5  32  0  32. 5  33.0  33.5  33.0  33.5  35.0 


SRLIIJIT. 

TEAR  NO.  1 
VARIABLE  MIXED  LAVER 
PIDJEx  CASE  3 

Figure  55.  Salinity  profile  from  May  10,  1975,  camp  Blue  Fox.  This  profile  is  used  to  initialize  the  VML  model. 
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Figure  56.  Mixed-layer  temperature  for  a  1-year  simulation  of  the  VML  mode l  forced  with  AIDJEX  forcing 
data.  Simulation  runs  from  May  I  to  April  29. 
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Figure  60.  Mixed-layer  salinity  measured  at  camp  Blue  Fox  from  May  10,  1975,  to  April  20,  1976. 
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Figure  61.  Mixed-layer  salinity  measured  at  camp  Snowbird  from  May  10,  1975,  to  April  20,  1976. 
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Figure  62.  Mixed-layer  depth  computed  from  a  1-year  simulation  of  the  VML  model  forced  with  AIDJEX 
forcing  data.  Simulation  runs  from  May  1  to  April  29. 
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Figure  63.  Mixed-layer  depth  measured  at  camp  Blue  Fox  from  May  10,  1975,  to  April  20,  1976. 
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Figure  66.  Comparison  of  temperature  profiles  from  Case  3  and  camp  Blue  Fox  from  September  !5  to 
January  I. 


87 


TEMPERATURE  PROEILE  pL0TS  EOR  MODEL  RUN  ONE  BLUEFOx  DRTR 
i ERR  NO.  1 


TENPEMTuPE  PROFILE  PLOTS  PER  NOlEl.  PUN  FIND  BLUEFOX  DflTR 
. ERR  NO.  1 


DEPTH  l METERS) 


100  90  80  70  60  50  10  30  20  10  0 


DCPTH  IMETCRSI 

100  90  80  70  60  50  10  30  20  10  0 


if* 

V /> 


Figure  69.  Comparison  of  salinity  profiles  from  Case  3  and  camp  Blue  Fox  from  September  15  to  January  1 . 
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Figure  71.  Comparison  of  temperature  profiles  from  Case  3  and  camp  Snowbird  from  May  IS  to  September  1. 
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Figure  72.  Comparison  of  temperature  profiles  from  Case  3  and  camp  Snowbird  from  September  14  to 
January  1. 
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Figure  76.  Comparison  of  salinity  profiles  from  Case  3  and  camp  Snowbird  from  January  IS  to  April  15. 
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Figure  80.  Comparison  of  mixed-layer  salinity  for  Case  3N  and  Case  30.  Open  water  occurred  in  both 
cases.  Case  3N  is  the  dashed  curve,  and  Case  30  is  the  solid  curve. 
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Figure  HI.  Comparison  of  ice  thickness  for  Case  3N  and  Case  30.  Open  water  occurred  in  both  cases.  Case 
3N  is  the  dashed  curve,  and  Case  30  is  the  solid  curve. 
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Figure  83.  Comparison  of  temperature  profiles  from  Case  30  and  camp  Blue  Fox  from  May  15  to 
September  l. 


104 
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Figure  84.  Comparison  of  temperature  profiles  from  Case  30  and  camp  Blue  Fox  from  September  15  to 
January  I. 
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Figure  86.  Comparison  of  salinity  profiles  from  Case  30  and  camp  Blue  Fox  from  May  15  to  September  I. 
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Figure  87.  Comparison  of  salinity  profiles  from  Case  30  and  camp  Blue  Fox  from  September  IS  to  January  1. 
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Figure  88.  Comparison  of  salinity  profiles  from  Case  30  and  camp  Blue  Fox  from  January  IS  to  April  IS. 
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SALINITY  PROFILE  PLOTS  FOR  MODEL  RUN  AND  BLUITOa  OATH 

YEAR  NO.  1 


Appendix  A:  Ocean  Model  Equations 


a.  Ice  velocity  equation 

The  ice  momentum  equations  (Eqs.  8  and  9)  can  be  solved  by  first  substituting  Eqs.  (25)  and 
(26)  for  the  water  stress  components.  The  equation  for  the  x-component  of  the  ice  velocity  then 
becomes 


3u .  , .  .  T 

5-1  -  f(v  -v  )  +  -ax 
ot  1  g  in 


+  5.5E-3 


(ug-ui) 


(v  -v. ) 

8  1 


((u^-u^cos©  - 


(vg-v^)sin©  ) 


(Al) 


where  S  =  25°. 

The  equation  for  the  y-component  becomes 

3v . 


^i_  -  -f(u.-u  )  *  ^ay  *  5.5E-3  1  1  [(u  -u.)'  +  (v  -v.)3] 

—  1  g  m  1-pihi-1  L  8  1  g  1'  J 

((Vg-vi)cose  *  (Ug-u.)sin©)  .  (A2) 

Next  let  C  =  5.5E-3  fp  / D  h  1  -J  [(u  -u  )’  +  (v  -v  )3]  and  rearrange  the  terms  in 
Eq.  (Al)  to  get  l  0  i  iJ  L  g  i  g  J 


_ i  *  (f+C  sin©) v  -  (C  cos©)u.-  fv  +  ^ax  +  C  u  cos© 

3t  1  1  S  m  g 

-  C  vgsin©  ,  (A3) 

and  in  Eq.  (A2)  to  get 


®vi  s  -  (f*C  sin©)u^-  (C  cos©)v.*  fu  +  J*ay  C  v  cos© 

3t 

+  C  UgSin©.  (A4) 

To  simplify  the  algebraics  let 

F  *  -fv  ♦  -^ax  *  C  u  cos©  -  C  v  sin©  , 
x  g  m  g  g  (A5) 


and 


F  « 

y 


C  v  cos© 

g 


C  UgSiip©  , 


f'  *  f  *  C  sin©  . 


(A6) 

(A") 


1 1 1 


Next  substituting  Eqs.  (A5),  (A6),  and  (A7)  into  Eqs.  (A3)  and  (A4)  yields 


3ui  «  f'v.  -  (C  cos©)  u.  +  F  , 

- —  1 

3t 


and 


3v . 


3 1 


«  -  f'u^  -  (C  cos0)  +  Fy  . 


The  finite  difference  form  of  the  above  equations  is  as  follows 


n+1  n 


n+1 


V  Ui  =  f'v.n+i  -  (C  cos0)u<n*1  +  Fvn  , 

— 5? -  1  1 


and 


(A8) 


(A9) 


(A10) 


n+1  n  ,  ,  _  n 

v.  -  v.  ,,  n+1  n+1  +  F 

_i _ i_  =  -f'u.  -  (C  cos0)v.  y 

At  1  1 


(All) 


Next  Eq.  (A10)  and  Eq.  (A  1 1)  must  be  solved  for  the  future  value  in  terms  of  the  previous  value, 
so  first  solve  for  u ln~l  and  for  The  above  equations  then  become 


and 


i  n  .  ,,  n+1  _  n 

n+1  u.  ♦  At  f  v.  +  At  F 
u  .  =  _ l _ l _ x 

1  1  +  C  At  cos0 


n+1  v.n+1-  At  f'  u.n+1+  At  F  n 

vi  =  __i _ i _ 1 


(A  12) 


1  +  C  At  cos© 

Finally,  substitute  Eq.  (A13)  into  Eq.  (A12)  and  again  solve  for  u,"'1  to  get 

n+1  (1+C  At  cos0)  u."+  (1+C  At  cos©)  At  F  n+  At  f'  v"+  At*f'  F  n 

u .  =  l  _  _ *  i  j 

1  2 

(1+C  At  cos 0)  +  At  f' 


(A13) 


(A14) 


The  y-component  of  the  ice  velocity  can  now  be  found  by  substituting  the  value  obtained  for  u l"*/ 
into  Eq.  (A13). 


b.  Water  velocity 

The  momentum  equations  for  the  ocean  (Eqs.  12  and  13)  are  solved  using  a  time-splitting  tech¬ 
nique.  First,  since  v  is  small,  it  is  combined  with  the  KM  term.  Next,  the  x-comoonent  of  the  veloci¬ 
ty  is  written  as  follows 


n  +  1 

u  -  u 
u _ w 

At 


n 


f  ♦v 


n  +  1  3  [km  [ 3u  nl  1 

/  +  £(1~*>Vv  *  >  J  ‘ 


f  v* 


(A15) 


where  $  can  range  from  0  to  1.  In  most  simulations,  f  =  1  is  used. 
Next  solving  Eq.  (A15)  for  u,'1*'  yields 


u 


n+1 


v 


uvn+  Atf^H*1*  Atf(l-*)v  H*1* 


At  f vn . 

g 


(A16) 


112 


Next,  an  intermediate  value  uw*  is  assumed  such  that 

it 


^  ^  • 


(A17) 


This  equation  can  be  solved  for  uw *  by  using  a  tridiagonal  solver.  Using  this  value  of  tr„  *, 
Eq.  (A  16)  can  be  rewritten  as 


n+1  *  , 

u  -u  , .  n+1  ...  ..  n  .  n 

_v _ v  «  f*v  +  f(l-#)v  -  fv 

At  8 


Then  as 

u  =  u  +  Atf*v  n+*  +  Atf(l-*)v  n-Atfv  . 

V  V  V  V  g 

In  a  similar  fashion,  the  y-component  of  the  velocity  can  be  written  as 

v  n*^-  v  ..  n+1  n  ,  n 

v  v  -  f*u  -  f ( 1 — ♦) u  +  fu_  . 

V  V  g 


At 


then  as 


(A18) 


(A  19) 


(A20) 


v  n+l=  v  -  Atf*u  n**-  At f ( 1  - d>) u  n  +Atfu  n  . 
v  v  v  v  g 


(A21) 


Finally,  substituting  Eq.  (A21)  into  Eq.  (A19)  and  solving  for  u ,nW  yields 
n-1  u  *♦  Atf *v  *-  (Atf)  2(l-«f)$u  "-Atf (l-4>)v  n+  (Atf)%u  n-  Atfv  n 

t  *  »  _  V*  V  V  g  g 


£ 


1  *  (Atf*)'  (A22) 

The  value  for  vw"*;  can  be  obtained  by  substituting  the  value  computed  for  uH"*;  into  Eq.  (A21). 


c.  Parameterization  of  turbulent  eddy  fluxes 

The  turbulent  eddy  fluxes,  KM  and  KH,  are  given  by  the  following  equations 
KH  ‘  lqSH  ’ 

and 

Km  *  IqSM. 


(A23) 


(A24) 


where  /  is  the  turbulence  length  scale,  q  is  the  square  root  of  twice  the  turbulent  kinetic  energy, 
and  SH  and  SM  are  stability  functions,  which  are  a  function  of  the  gradient  Richardson  number 
Ri,  where 


£  i£ 


Here.  g  is  the  acceleration  of  gravity  and  q  is  the  mean-field  density. 


11? 


The  quantity  q  is  calculated  from  a  form  of  the  turbulent  kinetic  energy  equation  that  expresses 
a  local  balance  of  shear  production,  buoyancy  production,  and  viscous  dissipation  of  turbulent 
kinetic  energy.  The  equation  is 

The  turbulence  length  scale  is  calculated  from  the  ratio 
turbulence  field  (Mellor  and  Durbin,  1975).  Thus 


151  =  0  ' 

(A26) 

of  the  first  to  the  zeroth  moment  of  the 


1 


0.1  | z |  q  dz 

J_mq  dz 


(A27) 


These  equations,  along  with  Eqs.  (12)  through  (15),  close  the  turbulence  parameterization. 


Appendix  B:  List  of  Symbols 


Symbols  Definition 
c  Specific  heat  of  seawater 

Coriolis  parameter 
Sum  of  the  atmospheric  heat  fluxes 
FB  Heat  flux  from  the  ocean  to  the  ice  or  atmosphere 

Fbr  Heat  stored  in  the  brine  pocket  reservoir 

F,  Latent  heat  flux 

Fl  Incoming  longwave  radiation 

FR  Incoming  shortwave  radiation 

FRH  Heat  from  solar  radiation  penetrating  through  the  ice  into  the  water  column 

Fs  Conductive  heat  flux  through  the  snow/ice  system 

Downward  flux  of  solar  radiation 
F.  Sensible  heat  flux 

Hdrtp  Heat  flux  from  the  deep  ocean 

hi  Ice  thickness 

h  mm  Minimum  ice  thickness  for  which  the  brine  pocket  simulation  is  calculated 

Maximum  heat  allowed  to  accumulate  in  the  brine  pockets  heat  reservoir 
hs  Snow  thickness 

10  Percentage  of  solar  radiation  that  penetrates  into  the  ice 

A',.y  Eddy  diffusion  coefficient  for  heat 

kt  Heat  conductivity  of  ice 

Av,  Eddy  diffusion  coefficient  for  momentum 

kt  Heat  conductivity  of  snow 

m  Mass  of  ice 

qb  Volumetric  heat  of  fusion  for  ice  at  the  bottom  of  the  ice 

q ,  Volumetric  heat  of  fusion  for  ice  at  the  top  of  the  ice 

qs  Volumetric  heat  of  fusion  for  snow 

S  Ocean  salinity 

S0  Reference  salinity  for  the  ocean 

t  Time 

T  Ocean  temperature 

TB  Temperature  at  the  bottom  of  the  ice 
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Ts  Temperature  at  the  surface  of  the  snow/ice  system  or  the  ocean  if  open 

water  occurs 

ua  x-component  of  the  wind  velocity 

ug  x-component  of  the  geostrophic  ocean  current 

w,  x-component  of  the  ice  velocity 

uw  x-component  of  the  ocean  current 

va  y-component  of  the  wind  velocity 

vg  y-component  of  the  geostrophic  ocean  current 

v,  y-component  of  the  ice  velocity 

vH  y-component  of  the  ocean  current 

w  z-component  of  the  ocean  current 

z  Vertical  coordinate,  positive  downward  from  sea  surface 

(  )  Ensemble  mean  for  equations  (10)-(17) 

f  ' )  Departure  from  above-defined  averages 

o,  Albedo  of  ice 

a  s  Albedo  of  snow 

Albedo  of  open  water 

A hb  Predicted  change  in  the  ice  thickness  at  the  bottom  of  the  ice-growth  or  decay 

A ht  Predicted  change  in  the  ice  thickness  at  the  surface  of  the  ice-melt  only 

A /i  s  Predicted  change  in  the  snow  thickness 

AS  Salinity  change 

A  r  Time  step 

ATS  Predicted  change  in  the  surface  temperature 

v  Background  vertical  eddy  diffusion  coefficient 

Qa  Density  of  air 

Density  of  ice 

qo  Reference  density  for  the  ocean 

qs  Density  of  snow 

a  Stefan-Boltzman  constant 

t  x-component  of  wind  stress 

r0>  y-component  of  wind  stress 

t(I  x-component  of  stress  between  the  ice  and  ocean 

thi  y-component  of  stress  between  the  ice  and  ocean 
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Appendix  C:  Corresponding  Calendar  Date  and  Model  Time  Step 


Dates  that  correspond  to  time  steps 
in  climatology  test  cases 


Date 

Time  Step 

Jan  1-10 

1-30 

11-20 

31-60 

21-30 

61-90 

Feb  1-10 

91-120 

11-20 

121-150 

21-30 

151-180 

Mar  1-10 

181-210 

11-20 

211-240 

21-30 

241-270 

Apr  1-10 

271-300 

11-20 

301-330 

21-30 

331-360 

May  1-10 

361-390 

11-20 

391-420 

21-30 

421-450 

Jun  1-10 

451-480 

11-20 

481-510 

21-30 

511-540 

Jul  1-10 

541-570 

11-20 

571-600 

21-30 

601-630 

Aug  1-10 

631-660 

11-20 

661-690 

21-30 

691-720 

Sep  1-10 

721-750 

11-20 

751-780 

21-30 

781-810 

Oct  1-10 

811-840 

11-20 

841-870 

21-30 

871-900 

Nov  1-10 

901-930 

11-20 

931-960 

21-30 

961-990 

Dec  1-10 

991-1020 

11-20 

1021-1050 

21-30 

1051-1080 

Dates  that  correspond  to  time  steps 
in  AIDJEX  test  cases 


Date 

Time  Step 

May  1-10 

1-40 

11-20 

41-80 

21-31 

81-124 

Jun  1-10 

125-164 

11-20 

165-204 

21-30 

205-244 

Jul  1-10 

245-284 

11-20 

285-324 

21-31 

325-368 

Aug  1-10 

369-408 

11-20 

409-444 

21-31 

445-492 

Sep  1-10 

493-532 

11-20 

533-572 

21-30 

573-612 

Oct  1-10 

613-652 

11-20 

653-692 

21-31 

693-736 

Nov  1-10 

737-776 

11-20 

777-816 

21-30 

817-856 

Dec  1-10 

857-896 

11-20 

897-936 

21-31 

937-980 

Jan  1-10 

981-1020 

11-20 

1021-1060 

21-31 

1061-1104 

Feb  1-10 

1105-1144 

11-20 

1145-1184 

21-29 

1185-1220 

Mar  1-10 

1221-1260 

11-20 

1261-1300 

21-31 

1301-1344 

Apr  1-10 

1345-1384 

11-20 

1385-1424 

21-29 

1425-1460 
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